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Abstract. Exponential Levy processes can be used to model the evolution 
of various financial variables such as FX rates, stock prices, etc. Considerable 
efforts have been devoted to pricing derivatives written on underliers governed 
by such processes, and the corresponding implied volatility surfaces have been 
analyzed in some detail. In the non-asymptotic regimes, option prices are de- 
scribed by the Lcwis-Lipton formula which allows one to represent them as 
Fourier integrals; the prices can be trivially expressed in terms of their implied 
volatility. Recently, attempts at calculating the asymptotic limits of the im- 
plied volatility have yielded several expressions for the short-time, long-time, 
and wing asymptotics. In order to study the volatility surface in required de- 
tail, in this paper we use the FX conventions and describe the implied volatility 
as a function of the Black-Scholes delta. Surprisingly, this convention is closely 
related to the resolution of singularities frequently used in algebraic geometry. 
In this framework, we survey the literature, reformulate some known facts re- 
garding the asymptotic behavior of the implied volatility, and present several 
new results. We emphasize the role of fractional differentiation in studying the 
tempered stable exponential Levy processes and derive novel numerical meth- 
ods based on judicial finite-difference approximations for fractional derivatives. 
We also briefly demonstrate how to extend our results in order to study im- 
portant cases of local and stochastic volatility models, whose close relation 
to the Levy process based models is particularly clear when the Lewis-Lipton 
formula is used. Our main conclusion is that studying asymptotic properties 
of the implied volatility, while theoretically exciting, is not always practically 
useful because the domain of validity of many asymptotic expressions is small. 
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1. Introduction 

In the classical Black-Scholes-Mcrton (BSM) European option pricing model (see 
[22] and [93]), asset processes are assumed to be strictly diffusive in nature and char- 
acterized by a single (log-normal) volatility a. In practice, no actual option market 
conforms with this framework, so to make the BSM formula work, practitioners are 
forced to make the volatility argument in this formula depend on option maturity 
(t) and strike (K). Indeed, it is common practice in virtually all option markets to 
maintain a strike- and maturity-dependent implied volatility surface, ai mp (t;T,K), 
such that a call option on an asset F paying (F(T) — K) + at expiration time T > t 
has a time t undiscounted price C(t; r, K) given by 

(1-1) ^|^ = $(rf + )- e fe $(d_), 

where k = \n(K/F(t)), r = T — t, $(•) is the cumulative Gaussian distribution 
function, and 

, , -k±h<rimv(t;T,K) 2 T 

(1.2) d± = 



Here and below, as usual, 

(x)* = ± max (±x, 0) . 

Note that our version of the BSM formula assumes that the asset price F is a risk- 
neutral martingale, an assumption that is easily relaxed or, in any case, justified if 
we consider F a forward process. At time t, the i-observable function <Ji mp (t; r, K) 
can be implied (with assistance of interpolation and extrapolation techniques) from 
quoted call option prices at multiple maturities and strikes. For later use, notice 
that the term 

(1.3) A = *(d+), 

is known as the option's (forward) delta. 

Below we often use the time value of a call option defined as follows 

(1.4) SC (t; r, K) = C(t: r, K) - (F (t) - K) + , 
or, equivalently, 

(1.5) 6C %f ] = - ^(d-) -{l-e k ) + . 

For future reference, it is convenient to introduce the following non-dimensional 
function C BS (v, k) of the annualized variance v and log-strike k: 

(1.6) C B5 (v,fc) = $(d + )-e fe $(d_), 
where < v < oo, — oo < k < oo, and, 

d± = 7 J-. 

In the limiting cases we have 

C ss (0,fc) = (l-e fe ) + , C BS (oo,fc) = l. 
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It is clear that 

C(t;r,K) BS f (f . TK ,2 T ,(K\ 
In most real markets, the market implied volatility function (Jimp 

(t; t, K) differs 

very significantly from a constant, and may have considerable slope and convexity 
as a function of K for both very small and very large values of r. To understand 
and explain this phenomenon, several alternative models have been proposed in the 
literature, see, e.g., [94], [39], [59], [14], [4], [76] among others. Broadly speaking, 
the available approaches can be categorized as follows: 

• Local volatility (LV) models, where the instantaneous volatility a of F is a 
deterministic function of time and F. 

• Stochastic volatility (SV) models, where a is a random variable, possibly 
correlated with F. 

• Jump models, where the process for F is assumed to be a purely discontin- 
uous jump-process. 

• Universal volatility (UV) models, where local volatility, stochastic volatility, 
and jump processes are combined. 

Full-blown UV models are rarely used in practice, and instead markets tend 
to converge around a simpler model that ultimately becomes a de facto market 
standard. For instance, LV models are popular in the field of equity derivatives, and 
jump-free combinations of LV and SV (known as LSV, or local-stochastic volatility, 
models) are dominant in the FX options arena. Such usage of simplified models 
is, however, rarely motivated by empirical facts, but instead are done for practical 
reasons in recognition of the fact that models that combine LV, SV, and jumps 
are often highly complicated to calibrate and implement. One particularly thorny 
issue is the question of how precisely to mix LV, SV, and jumps, a problem that 
is made especially vexing by the fact that even a simple model like LV can, on its 
own, match essentially any arbitrage- free option volatility surface (see, e.g., [39]). 

In UV models, there are several potential strategies for attacking the "mixing" 
problem. As LV, SV, and jumps give rise to different dynamic evolution of the 
volatility surface over time, one idea is to embed observations of the volatility dy- 
namics into the model calibration problem (see, e.g., [6] for a discussion). A closely 
related approach is to incorporate market-observable exotic option prices that are 
sensitive to the evolution of volatility, e.g. variance options, barrier options, and 
similar; this approach is commonly used in FX markets to mix LV and SV into 
an LSV model. Yet another idea is to examine various extremes of the volatility 
surface (short maturity, long maturity, large strikes, small strikes) and attempt to 
understand which type of model feature (SV/LV/jumps) will control the asymp- 
totic model behavior most strongly and most "naturally" . This information, in 
turn, could then be used to inspire the model building process. To give an exam- 
ple, consider that the convexity of <7j TOp (t;T, K) around K = F(t) can often be 
observed to be substantially higher for small option maturities than for larger ones. 
A pure LV model would model this by letting the convexity (in F) of the local 
volatility function decay rapidly as a function of time. Such a model would, how- 
ever, be highly non-stationary, an undesirable model feature that can be avoided 
by introducing jumps or (mean-reverting) stochastic volatility. 
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The idea of considering volatility function asymptotics is particularly attractive 
from an analytical perspective, as one can often work out simple closed-form ex- 
pressions for these asymptotics, even for complex models. For SV, LV, and LSV 
models the procedure for generating such asymptotics is well-understood (see e.g. 
[58] for a survey of recent results), and the range of validity for the resulting expan- 
sions has been examined closely, generally with fairly satisfying results. Recently, 
similar attempts have been undertaken for jump processes, especially for processes 
in the so-called exponential Levy class. More precisely, these models write 

(1.7) F(t) =F(0)e* (t) , 

where X(t) is a Levy process. In this setup, particular emphasis has been put 
on the small-time asymptotic of C(t;r,K) for r — > and K fixed at some level 
different from F(t), see for instance [70], [105], [46], [48] among many others. A 
typical result in this area of research is that 

(1.8) 5C(t;r,K)~T, K £ F(t), 

where ~ indicates the leading order term as t — > 0. For this relation to be true, 
the corresponding implied volatility must explode in the short-time limit. This line 
of research is closely related in studying small time asymptotics for the densities 
of Levy processes, see, e.g., [67], [98], [107], [44], [45], [84], [18]. Small-time results 
for at-the-money (ATM) options where K = F(t) are scarcer, but in [30] Tanaka's 
formula is used to demonstrate that 

(1-9) C(t;r,F(t))~{^ -J 

when X is a mixed jump-diffusion consisting of a Brownian motion with volatility a 
and a finite variation 1 Levy jump process. [116] further elaborates on this result in 
implied volatility space, and demonstrates that for a finite variation Levy process 
(necessarily without a diffusion component) 

(1.10) <i mp (t; r, F(t)) ~ V^max ( / (e x - 1)+ v (dx) , [ (1 - e x ) + v (dx)\ , 

Ur Jr J 

where v{-) is the Levy measure of X. In the presence of a diffusion component with 
constant volatility a, [116] shows that if 



/ 

Jr 



x v{dx) < oo, 

Jr 
then 

Pimp {t;T,F(t)) ~ a 

for t — »■ 0. Other relevant papers about short-time option prices in models with 
jumps include [2], [91], and [96]. 

As for long-maturity and extreme strike asymptotics, the literature is generally 
quite sparse. Under mild regularity conditions, [117], [103] use large deviation 
principles to demonstrate that the implied volatility surface will always flatten out 
for large enough maturities and finite strikes, irrespective of the underlying process 
assumptions, which agrees with earlier observations by [11]. For exponentiated Levy 
processes, the authors list explicit formulas tying the long-time volatility asymptote 
to moments of F. Extreme strike asymptotics for the specific case of Levy processes 

■'■See Section 3 for a definition of various attributes of Levy processes. 
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follow from the moment formulas of [69], [15], [16] and [55]; see [116] for a brief 
discussion. 

Overall, while some progress has been made in recent years, the available asymp- 
totic results for jump processes are limited when compared to the literature for LV 
and SV diffusion processes. In particular, little is known about the range of validity 
of jump process asymptotics, in part due to the difficulty of establishing accurate 
numerical option prices for Levy processes, especially in the short-time limit. In 
this paper, we add to the literature in several ways. First and foremost, we com- 
pletely characterize the r and K volatility surface asymptotes for the important 
class of (exponentiated) tempered stable Levy processes, with and without a dif- 
fusive component. In addition, we also list a number of asymptotic expansions for 
the Merton jump-diffusion and the NIG (Normal Inverse Gaussian) process. Note 
that our results cover not only the level of the volatility surface, but also its slope 
and convexity - properties that are highly relevant when calibrating a model to ob- 
served option prices 2 . We mainly work in Fourier space using the Lewis-Lipton call 
option representation [72] , [76] , and our primary tools for asymptotic expansion are 
the dominated convergence theorem (for the short-time limits) and classical com- 
plex variable techniques, such as the saddlepoint theorem (for long-dated options) 
and high-frequency Fourier asymptotics (for large and small strike limits). For 
testing and illustration purposes, we present several closed-form option prices for 
special cases of the tempered stable Levy process, and also list some results for the 
computation of Green's function. Further, we draw attention to the applicability 
of fractional differentiation in characterizing the PIDE (partial integro-differential 
equation) that governs option prices in the tempered stable model class; these re- 
sults allows one to draw on modern numerical algorithms when pricing both vanilla 
and exotic options. As an example, we develop an operator splitting method which 
is O (N 2 ) complex and second order accurate in time and space. The idea is to 
split jumps into positive and negative and treat them separately on each step; the 
problem then boils down to the inversion of Hessenberg matrices which may be ac- 
complished via the generalized Thomas algorithm in O (iV 2 ) operations. A regular 
diffusion can be added as an additional iteration, as needed. 

To test and motivate our asymptotic results, our primary focus in this paper is 
on FX options markets. These markets are of interest to us for two reasons: i) very 
short-dated options actively trade in this market (it is not uncommon for dealers to 
quote on options with a maturity being just a few hours); and ii) volatility quoting 
conventions in these markets use the option delta A (see (1-3)), rather than the 
strike K. At first glance the delta quotation standard may appear to be a nuisance, 
since additional translation from the delta space into the strike space is needed, but 
on second thought it becomes apparent that this construct has deep mathematical 
roots. Specifically, quoting volatility as a function of delta is closely related to the 
so-called resolution of singularities frequently used in algebraic geometry and other 
mathematical disciplines (see [9]). Our tests disprove a number of "urban myths" 
about Levy process asymptotics, especially regarding the range of applicability of 
short- and long-time asymptotics. 

Our paper is organized as follows. Section 2 provides a very brief introduction to 
FX volatility quotation standards, and introduce the concepts of risk reversals and 



Indeed, in FX markets it is standard to effectively quote directly on smile slope and convexity, 
see Section 2. 



EXPONENTIAL LEVY PROCESSES 



7 



straddles (also known as butterflies). We examine some representative market data 
for implied volatility, and highlight how the short-dated asymptotics are unnatural 
in a diffusion setting. Section 3 gives an overview of Levy processes, with an 
emphasis on the tempered stable class and the fractional derivative representation 
of the corresponding option pricing equation. It also discusses both conventional 
and novel numerical methods for the tempered stable class which capitalize on their 
interpretation as fractional derivatives. Section 4 discusses some representative 
examples. Section 5 introduces the Lewis-Lipton formula which is the main working 
tool for establishing all our asymptotic results. It also discusses calibration of 
the corresponding stochastic processes to market data introduced in Section 2. 
Section 6 discusses the asymptotic behavior of the BSM formula which is used 
subsequently in order to analyze the behavior of implied volatility. Sections 7, 8 
study long-time and short-time asymptotics, respectively, and Section 9 deals with 
wing asymptotics. Besides deriving theoretical results, Sections 7, 8, 9 also contain 
a series of numerical tests aimed at establishing their practical relevance. Finally, 
conclusions are drawn in Section 10. Appendices contain some of the more elaborate 
proofs and other useful material. 

2. Background on FX Market 

The FX options market is one of the largest over-the-counter options markets 
in the world, yet its conventions are quite idiosyncratic and differ markedly from 
those used in other derivatives markets (e.g., interest rates and equities). Moreover, 
almost every concept of importance can be interpreted differently, often depending 
on the currency pair in question, which makes systematic analysis and comparison of 
FX options particularly difficult. There are both historical and financial reasons for 
the existence of FX quotation styles, the most obvious being that FX transactions, 
unlike many other financial transactions, are inherently symmetric in nature, in the 
sense that units of currency are exchanged into units of currency (see, e.g.. [75]). 

In FX markets, it is standard practice to represent the volatility smile in terms 
of the option delta, rather than in terms of the option strike. We may define 

(2.1) S(t;r,A) = a imp {t;T,K) 

where A is the delta defined in (1.3). The map between A and K is monotonic, so 
(2.1) is always well-defined. We should note that several other definitions of delta 
than (1.3) exist in the FX market, not all of which are monotonic in strike 3 ; for our 
purposes, we ignore this complication and just refer to [102] and [33], among many 
other sources, for additional information on various delta definitions. 

In the FX options markets, the function E(t; r, A) is normally liquidly quoted 
at only three different levels of delta: 0.25, 0.50, and 0.75. Somewhat confusingly, 
only X(i;r, 0.5) is directly quoted (the at-the-money volatility 4 <jatm, where K = 
F(t)), whereas £(i;r, 0.25) and E(i;r, 0.75) must be constructed from quotes for 



For instance, the so-called premium- adjusted (forward) delta &(d-)K/F is not mononic in 

K. 

4 In real FX options markets, the ATM strike may actually deviate slightly from K = F, as 
additional conventions govern the choice of the at-the-money strike. For instance, a common 
alternative is to use the delta-neutral strike K = F cxp(crj mp T /2) , which is the strike level where 
the absolute magnitude of put and call deltas are identical. 
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risk-reversals (RRs) and butterflies (BFs) (also known as strangles). The relevant 
formulas 5 are 



BF(t;r) 



S(i;r,0.25), 
(£(i; r, 0.75) + r, 0.25) 



2& atm^t)) , 



oatm^t) = £(t;r, 0.5), 
RR{Ut) = E(i;r,0.75) 
1 
2 

which trivially allows us to construct £(i;r,0.75) and £(i;T,0.25) from knowledge 
of (J atm , -R-R, and Bf . It is clear that the RRs and BFs are closely related to the 
slope and convexity, respectively, of r, A) around A = 0.5. Specifically, we can 
write (omitting arguments) 

i as 1 9 2 S 

(2 ' 2) RRRi 2dA ; BFKS 32d&' 

where the derivatives are evaluated around A = 0.5. 

In Table 1 below, we show some sample market quotes for a atm {t), RR(t), 
and BF (t) for the USD/JPY currency pair. We highlight that the figure suggests 
the existence of non-zero finite limits for all three quotes a atm, RR, and BF as 
r approaches zero. From (2.2), this essentially translates into non-zero finite limits 
for both <9£/<9A and <9 2 £/<9A 2 around the ATM point. 





Id 


lw 


2w 


lm 


2m 


3m 


6m 


iy 


2y 


3y 


4y 


5y 


ATM 


9.87 


12.15 


11.27 


10.35 


10.34 


10.52 


11.04 


11.80 


12.90 


13.80 


14.30 


15.05 


RR25 


-1.25 


-1.00 


-0.46 


-0.30 


-0.16 


-0.10 


-0.05 


-0.02 


-0.35 


-0.78 


-1.20 


-1.48 


BF25 


0.30 


0.30 


0.30 


0.30 


0.32 


0.33 


0.46 


0.59 


0.60 


0.58 


0.53 


0.46 


RR10 


-2.22 


-1.77 


-0.79 


-0.49 


-0.22 


-0.09 


0.07 


0.18 


-0.50 


-1.41 


-2.30 


-2.86 


BF10 


1.11 


1.08 


1.03 


1.02 


1.08 


1.12 


1.54 


1.98 


2.11 


2.22 


2.28 


2.29 



Table 1. The behavior of <tatm{t), RR(t), and BF (t) (ex- 
pressed in per cent) for USD/JPY on March 30th, 2012. RR25 
and RR10 are the risk-reversals at A = 0.25 and A = 0.10, respec- 
tively (and similar for BF25, BFIQ). 



To understand the implications of Table 1 for the volatility smile in strike space, 
recall our definition of log-moneyness k = \n(K/F(t)) and notice that 

ln „x <9X d<Ti mv r- — d 2 T, I ' d 2 Oi mri 1 d(Ji mD \ „ 

where the derivatives are taken around A = 0.5 or, equivalently, k — 0. To match 
non-zero limits of 9£/<9A and <9 2 £/9A 2 for t — > 0, we evidently need both the 
smile skew d<Ji mp /dk and the smile convexity d 2 a i mp / dk 2 to approach infinity at 
k = for small r, at rates of t~ x / 2 and r _1 , respectively. Such a requirement, 
however, would rule out that FX dynamics are driven by a pure diffusion process 
(such as LSV models), since it is known (see Section 8.3) that such processes always 

5 We are omitting some complications in the definition of risk reversals and butterflies. Besides 
being imprecise about the correct definitions of delta and the ATM strike, we have chosen to define 
the butterflies as so-called smile stranges, rather than the more common market strangles. The 
latter definition can be found in [33] and, unfortunately, does not allow one to uniquely extract 
E. 
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result in a finite limit of dai mp /dk at k = 0. Motivated by this observation, we 
proceed below to introduce the class of Levy processes. 

The data in Table 1 can be converted (with the aid of interpolation and extrap- 
olation) into the function E (r, A). The result is shown in Figure 1. It is clear that 
the FX options exhibit strong smile for short and medium maturities. 




Maturity 



Figure 1. USDJPY implied volatility surface. The quotes are for 
March 30th, 2012. 



3. Background on Exponential Levy Processes 

3.1. Basic Setup. In this section we consider exponential Levy processes (ELPs). 
General properties of these processes are discussed in detail by numerous researchers, 
see, e.g., [21], [109], [8], among many others. The realization that these processes 
have important applications to finance can be traced back to [83] (and, in fact, to 
much earlier work); these applications are discussed in many books and papers, see, 
e.g., [41], [17], [101], [13], [24], [111], [7], [34], to mention just a few. 

Let X(t) be a Levy process, i.e. a cadlag process with stationary and independent 
increments, satisfying X(0) = 0. It is known that every Levy process is characterized 
by a triplet (7,er, is), where 7 and a axe constants, and v is a (possibly infinite) 
Radon measure, known as the Levy measure. The Levy measure must always satisfy 

Jr 

To characterize the infinitesimal generator of a Levy process, let E be the expecta- 
tion operator, and define 

V(t,x)=E(V T (X(T))\X(t)=x) 
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for some suitably regular function Vr(-). It can be shown that V solves a partial 
integro-differential equation (PIDE) of the form 

V t +iVx + \<r 2 V xx + [ (V(t,x + y)-V(t,x)-yl lvl < 1 V x (t,x)) v{dy)=0, 
1 Jr 

subject to the terminal condition V {T,x) — Vt (x). More conveniently, we write 

this PIDE as 

(3.2) 

V t +'yV x + la 2 (V xx -V x )+ [ (V(t, x + y) - V (t, x) - y\\ v \<iV x (t, x)) v(dy) = 0, 

1 JR 

where 7 = 7 + cr 2 /2. Loosely speaking, (3.2) demonstrates that any Levy process 
X(t) is the sum of a deterministic drift term yb, a scaled drift-adjusted Brownian 
motion aW(t) — a 2 t/2, and a pure jump process. 

By choosing Vr {x) = exp(iux) and solving (3.2) through an affine ansatz 

V(t, x) = exp(tp(u)(T -t) + iux), 

one arrives at the famous Levy-Khinchine formula, 

u) = E = cxp (V>(u)t) , 

(3.3) ip(u) = Jin - -a 2 u(u + i) + [ (e mx - 1 - iuxl^^) v(dx), uGR, 



2 

where ip is the so-called Levy exponent. In practice, a Levy process may be specified 
either by its exponent ip or by its Levy measure v. It is frequently convenient to 
split ip(u) into two parts as follows 

(3.4) i>{u) = - l -o- 2 u{u + i)+il) {u), 

where ip (u) is the pure-jump and drift component of tp(u). 

Let F(t) > denote the time t price of a financial asset, and consider modeling 
its evolution as an exponential Levy process of the type (1.7), where we wish for 
F(t) to be a martingale in some pricing measure. For this, we impose that, for any 
t > 0, 



(3.5) E 



1. 



which requires that the first exponential moment of X(t) exists in the first place, 
i.e. that large positive jumps be suitably bounded: 

(3.6) / e x v{dx) < 00. 

J\x\>l 

Equivalently, we require that <j>(t, z) exists in the complex plane strip 

S = {ze C|Im(z) G [-1,0]}, 

and that 

4>(t, -i) = cxp (jt + t j (e x - 1 - zl| x |<i) v(dx)j = 1, 

in order to satisfy condition (3.5). This implies the fundamental martingale con- 
straint 

(3.7) 7 = - / (e* - 1 - xh xl<1 ) u{dx). 
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Using (3.7) to eliminate 7 from (3.3), the martingale restriction on F allows us 
to write the Levy exponent in the form 

(3.8) V(u) = -\<? 2 <u + [ (e lux -l-iu{e x - 1)) v(dx). 

Similarly, we may write the PIDE (3.2) as 
(3-9) 

V t + \o 2 (V xx -V x )+ [ (V(t, x + y)-V (t, x) - ( e y - 1) V x (t, x)) v{dy) = 0. 
A Jr 

Whenever possible (see Section 3.2), it is often more convenient to use the fol- 
lowing forms of the PIDE (3.2), namely, 

1 f°° 
(3.10) V t +iV x + -a 2 (V xx -V x )+ (V(t,x + y)-V(t,x))u(dy)=0, 
1 J -00 

(3.11) 

1 f°° 
V t + i'V x + -a 2 (V xx -V x )+ (V (t,x + y)-V (t, x) - yV x (t, x)) v{dy) = 0, 

^ J -co 

where 

(3.12) 7' = - / (e x - 1) v{dx) = 7 - / xv(dx), 



JR J|a;|<l 

(3.13) 7" = - f (e x -l-x) v{dx) = 7 + f xv{dx). 

Jr J\x\>i 

3.2. Classification of Exponential Levy Processes. Depending on how "sin- 
gular" the Levy measure v is, we can define various sub-groups of Levy processes. 
Each group allow us to decompose equations (3.3) and (3.2) in slightly different 
ways. 

3.2.1. Finite Activity. First, if the Levy measure is finite (i.e., the jump component 
of the process has finite activity), the resulting process for X(t) is a combination 
of a Brownian motion and an ordinary compound Poisson jump-process. We may 
then replace v(dx) with 

(3.14) u(dy) = Xj(dy), A = / u{dx) < 00, 

Jr 

where j(dx) = v(dx)/\ is now a properly normed probability measure for the 
distribution of jump sizes in X, and A is the (Poisson) arrival intensity of jumps. 
In this case, 

i>(u) = j'iu - \o 2 u (u + i) + A ( / e iux j( dx) - 1 

where the martingale restriction requires that 7' satisfies (3.12). Notice that if we 
define a random variable J with density j(dx), then we can, in the finite activity 
case, interpret 

(3.15) ip{u) = i%u - ^a 2 u (u + i) + A {<t>j{u) - 1) , 
where 4>j{-) is the characteristic function of the jump size 6 J. 



^Specifically, if a jump of size J takes place at time t, X(t—) goes to X (t+) = X(t— ) + J and 
S(t-) goes to S(t+) = S(t-)e J . 
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As for the PIDE (3.2), for finite activity processes it simplifies in an analogous 
way to the simpler form (3.10): 

V t + iV x + \a 2 (V xx - V x ) + X f V(t, x + y) j(dy) - XV = 0. 
1 Jr 

3.2.2. Finite Variation. The jump component of a Levy process is said to have 
finite variation if 

(3.16) / \x\v{dx) < oo. 

J\x\<l 

Under this condition, truncation of the Levy exponent tp around the origin is not 
necessary, and we may write 

(3.17) ip{u) =iiu- \o 2 u{u + i)+ [ {e mx - l) v{dx), i = 7 - / xu{dx). 

2 Jr J\x\<i 

For the PIDE (3.2), we then get the simpler form (3.10). 

3.2.3. Finite First Moment. Finally, for the case where the first moment exists, 

(3.18) / {xll^yxvi dx) < 00, 

JR 

we may write the Levy exponent in purely compensated form: 

ip(u) = i'iu - \o 2 u (u + i) + f {e lux - 1 - iux) v(dx), 7" = 7 + / xv(dx). 
2 Jr J\x\>i 

The corresponding PIDE can then be written in the simpler form (3.11). 

Going forward we shall often omit primes on 7 and simply use 7 to denote the 
drift term of the PIDE, implicitly choosing the right one. 

4. Examples of Exponential Levy Processes 
4.1. Tempered Stable Processes. 

4.1.1. Definitions and Basic Facts. Establishing short-time ATM volatility smile 
asymptotics for the completely generic class of exponential Levy processes appears 
to be a difficult problem, so we narrow our focus to classes of processes important 
in applications. Of primary importance to us are processes characterized by Levy 
measures of the form 

(4.1) v{dx) = (^e-^Wo + ^T^ K - N ^<o) **, 

where we require 7 that n + > 1, k_ > 0, c + > 0, c_ > 0, and a < 2. The resulting 
class of processes 8 is known as tempered a-stable Levy processes (TSPs), see, e.g., 
[64], [89], [23], [24], [25], [34], [106]. Occasionally, a different parametrization of the 
Levy measure is used 

(4-2) v(dx) = - S ^M (J^-^Wo + ^- K - lX ^<o) *», 

7 The condition k+ > 1 is required to satisfy (3.6). In most literature, the condition is the less 
tight K-|- > 0. 

8 Only for a > does the class behave like a stable process, but we here allow for a < to 
include compound Poisson processes. 
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fl + > 0, i9_ > 0, and sec (n/2) /T(—l) = —2/ir. This parametrization is particularly 
useful when one analyzes changes occurring when a crosses unity. Asymmetry of 
the TSP is often characterized by the non-dimensional number (3, 

c + — c_ 



(4.3) 



c+ + c_ 

It is clear that TSPs are natural extensions of the classical a-stable Levy processes 
(SPs), where k± = 0, see, e.g., [122], [108], [97]. 

The overall behavior of the TS class is closely tied to the selection of the power 
a, as demonstrated in Table 2. 



a 


Activity 


Variation 


< 
(0,1) 
[1,2) 


Finite 
Infinite 
Infinite 


Finite 
Finite 
Infinite 



Table 2. The behavior of the TS Levy class as a function of a. 



Some important special cases of the TS class include: 

• The KoBoL (CGMY) model, where c+ = c_, see [64], [89], [23], [25]; 

• The exponential jump model, where a = —1, sec [65], [78]; 

• The Gamma process, where a = and either c_ = or c + = 0; 

• The Variance Gamma model, where a = 0, see [82]; 

• The Inverse Gaussian process, where a = 1/2 and either c_ = or c + = 0, 
see [12], [13]. 

For some of the special cases above, explicit formulas exist for the density of 
X(t) and for European call options. Section 4.1.4 lists such formulas for the Inverse 
Gaussian process, which we use for testing various results later on. 

When a^0 and a^l, the characteristic function for the TS Levy process can 
easily be shown to be, (see, e.g., [34],) 

(4.4) V(w) = r c * (( K * _ siu T ~ iu ( K s _ sl T + ( iu - *) O 

s=± 

— a s (k s — siu) a + ^iu + S, 

s=± 

where 

(4.5) a, = r(-a)c 8 = -sec(^)C C s = « - («. - *D , 
V s = 7 = a+C+ + a-C-, S = {a + rj + +a_r/_) . 

We have imposed the martingale condition (3.7) to express 7 as function of other 
parameters. In (4.4) the complex power functions 

(k_I_ — iu) a , (k,- + iu) a , 

are here (recall that a 7^ 0, 1, 2) multi-valued functions, and an appropriate branch 
cut is required. We need, as a minimum, that Tp{u) is regular for the strip S, i.e. 
when u = v! + iu", v! e R, u" e [—1,0]. In this strip, the two power functions 
evaluate to 

(k + + u" — iu') a , (k_ — u" + iu') a . 
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As the real part of the argument of the power function is strictly positive in <S, a 
branch cut in the left half-plane (e.g., the usual principal value branch cut for the 
logarithm) will therefore suffice. 

Equation (4.4) holds only for a ^ 1 and a^O (and a < 2, of course). When 
a = 0, we get 

<«» *'«'-£* K^)- ,og (;r^H 

Proceeding as above, we can easily show that the arguments to the log-functions in 
this expression are entirely in the right half-plane when Im« e [—1,0]; the principal 
value of the logarithm will suffice. Finally, for the case a = 1, we have 

ip{u) = ^ c s (k s - siu) log ^ - — ^ - (k s - s) log ^ Ks ^ iu. 

Again, we may interpret log as the principal value of the logarithm. 

By using (4.4) it is easy to show that the annualized standard deviation of a 
TSP has the form 

(4.7) stdev(X(l))=(r(2-a)(^ + -^ 

Equation (4.7) allows us to get an idea of the magnitude of the implied volatility 
of ATM options on TSPs. 

4.1.2. PIDEs and their Fractional Derivative Interpretation. Consider now the TS 
Levy class with an added Brownian motion with volatility a. For a e (0, 1) the 
PIDE (3.10) applies, and has the form 

(4.8) V t + iV x + -a 2 (V xx -V x ) + J2cs (V (x + sy) - V (x)) = 0. 

1 s=± J o y 

For a e (1, 2) the PIDE (3.11) can be used, 
(4.9) 

V t +i'V x +-o- 2 (V xx - V x ) + Y] c s / (V(x + sy) - V (x) - syV x (x)) y - 0. 

2 »=± y 

Interestingly, it is possible to rewrite both (4.8) and (4.9) in terms of so-called 
fractional derivatives (see, e.g., [95], [100] for a survey, and [122] for a connection 
to regular, non-tempered stable Levy processes). The development of fractional 
derivatives originated in the nineteenth century with Riemann, Liouville and Mar- 
chaud among others (see, e.g., [85]), and traditionally starts with the well-known 
Riemann-Liouville formula which allows one to reduce the calculation of an n-fold 
integral to the calculation of a single convolution integral. This formula can be 
extended in a natural way for a- fold integrals for any a > 0. Its extension for 
negative a, which can be viewed as a-fold differentiation, can be done in several 
different ways; we find that the so-called Caputo definition is the most convenient 
for our purposes. For integer values of a, the Caputo derivative coincides with the 
regular derivative of order a. For non-integer values of a, consider first a <G (0, 1), 
and define left (s = —1) and right (s = 1) fractional derivatives of order a as follows 
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®sV (x) = f (V (x + sy)-V (x)) , s = ±. 

Under some mild regularity assumptions we can perform integration by parts and 
write 

(-s)" -1 [°° dV (x + sy) dy 



»™ r<i-«)y„ * y 

For all non- integer values of a G (1, oo) we may then define 

(-s)"- 1 -^ f°° d^ +1 V (x + sy) dy 



V«V(x) = 



[I + [a\ - a) J 



r (1 + [a\ - a) J dyL«J+i ' 

where [a] is the floor function, i.e., the largest integer such that [a\ < a. Notice 
that in general D"V (x) is complex-valued even when V (x) is real-valued. 
We emphasize that with this definition, irrespective of s, 

®%e iux = (iu) a e iux , 

consistent with what one would expect from a generalization of a regular derivative. 
Lemma 4.1. For a e (0,2), a ^ 1, tfte P/L'Ss (^.5J, ^.Pj may be written 

(4.10) U t + jV x + \a 2 (V xx -V x ) + J2 a s e SK ° x V« (e^V) +5V = Q, 
where a s ,~f,5 are given by expression (4-5). 

Proof. See Appendix A. □ 

Corollary 4.2. In particular, for non-tempered stable processes with k± — 0, the 
corresponding PIDE has the form 

v t + \o 2 (v xx -v x ) + Y, (~ s ) a a s (®"V - V x ) = 0. 

1 s=± 

4.1.3. Numerical Methods Based on Fractional Derivatives. The advantage of re- 
stating the original pricing PIDEs as fractional differential equations is that we 
may lean on a large body of literature on numerical methods for such equations. 
These methods have been developed over the last twenty years due to the fact that 
fractional differential equations have numerous physical applications, especially to 
the so-called anomalous diffusions, see, e.g., [31], [112] among many others. In this 
section we present some extensions of these methods for our setting where both left 
(s = — 1) and right (s = +1) derivatives must be considered simultaneously. 

Before proceeding, let us remind the reader that an N x N matrix H such that 
hij = when j < i — 2 (hij = when j > i+2) is called an upper (lower) Hessenberg 
matrix. Such a matrix can be viewed as a generalization of a tri-diagonal matrix. 
An equation of the form 

(4.11) Hp = q, 

can be solved for the vector p via an appropriate extension of the Thomas algorithm 
for tri-diagonal matrices at the cost of O (N 2 ) operations. The advantages of using 
Hessenberg structure of the problem arc manyfold (see [114]), the most obvious 
being the ability to rely on highly parallclizable solvers. Hessenberg matrices nat- 
urally arise when one wants to solve pricing equations with fractional derivatives. 
When such derivatives are one-sided (or, equivalently, jumps are one-sided) we can 
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naturally represent them via Hessenberg matrices on a grid by virtue of an appro- 
priately discretized 55". When they are two-sided, we can split the problem into 
two and consider positive and negative jumps in turn at each step. Diffusion can 
be added as needed. The Grunwald-Letnikov formula is often used for discretizing 
fractional derivatives, see, e.g., [92], [38], [115]. However, it is more convenient to 
discretize 55" directly, see, e.g., [113] for a discretization scheme that guarantees 
that the resulting finite difference scheme is stable. In this discretization scheme, 
fractional derivatives 55" in (4.10) turn into Hessenberg matrices "H"; while the 
diffusion-advection term turns into a standard tri-diagonal matrix V. In order to 
take full advantage of the nature of the matrices Hg,T>, we solve a generic evo- 
lution equation of the form (4.10) by splitting a typical time step t m+1 — > < m , 

ym+1 yra ^ thm ^ ym+1 V* , V* ^ V** , V** -> V m as follows: 

I - l&tH%J V* = + l&tH%J V m+ \ 
I - ^AtH°^j V** = (j+ \ Am °^j V*> 

I - \&tv^j v m =(i+ \ AtT ^j v**- 

It is clear that at the first and second intermediate steps we have to solve a Hes- 
senberg system of equations, while at the third step we need to solve a tri-diagonal 
system of equations. We illustrate the scheme by numerically constructing the 
probability density function (PDF) of a typical TSP in two different ways, namely 
by solving the corresponding PIDEs and by numerically calculating its integral 
representation. Our results are shown in Figure 2, and it is clear that our method 
reproduces the corresponding PDF very accurately. 

There are several well-known Fourier-transform-based and approximation-based 
approaches to solving the corresponding pricing equation in the spatially homoge- 
neous case, see, e.g., [4], [10], [36], [32], [120], [60], [79], so, if this was the only 
case one is interested in, it would not be worthwhile adding yet another method 
to this set. However, to the best of our knowledge, none of these methods can 
be straightforwardly extended to the inhomogeneous case, nor to the case when 
barriers are present. Our method on the other hand can be extended to cover these 
cases almost verbatim. 

4.1.4. Maximally Skewed Stable Processes. If we set k + = k_ = 0, the Levy mea- 
sure (4.1) becomes that of a regular stable process. This type of process has limited 
applications in finance, as we only can satisfy the martingale restriction (3.5) if we 
mandate that c + = (/? = — 1, where j3 is given by (4.3)), i.e. when all jumps 
are negative. The resulting process is sometimes known as the maximally nega- 
tively skewed stable process, and has received some interest in the literature, see 
[29] . The characteristic exponent for such processes is known in closed form. Using 
parametrization (4.2) and assuming that the process is a martingale, ip(u) can be 
defined as follows 

[4.12) W-j |^nln(m), a = 1, 

in the complex strip S defined by Im(u) e [—1,0]. Despite the fact that jumps 
can only "go down" when (3 = — 1, the already present drift-correction around the 
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Figure 2. Numerical and analytical PDFs for a TSP. For com- 
parison, PDFs for the moment-matched Gaussian process is shown 
as well, see (4.7). The corresponding parameters are a = 1.50, 
a = 0.01%, c+ = 0.0069, c_ = 0.0063, k + = 1.9320, k_ = 0.4087. 
The choice of parameters is justified in Section 5.5. 



origin in the Levy-Khinchine theorem makes the resulting process a martingale for 
q = 1 , so no additional drift correction is necessary. 

Remark 4.1. There are three classical cases when the PDF for the stable Levy 
process can be calculated explicitly: i) a — 1/2,(5 — ±1 (Levy distribution); ii) 
a = l,/3 = (Cauchy distribution); Hi) a = 2, 1 < /3 < 1 (Gaussian distribution). 
Of these three cases, for option pricing purposes the only nontrivial case is a = 1/2 
and j3 = — 1, where 



V 2<d 

"( da; ) = ^7 -^H^-x<odx. 

2y/Tr(—x) 1 

Here, the distribution of X{t) is 

exp 

(4.13) P W*) e [x,x + dx]) = — 



A(l-x) 



^-l x<L dx, 



(l-x) 3/2 

where i = \/2~dt and condition (3.5) is clearly satisfied, so that exp(X(t)) is a 
martingale. Notice that in this case X(t) is bounded from above by l, which may 
not be particularly realistic. Still, for testing purposes it is worthwhile developing 
this case further, so Proposition J^.S develops an analytic call option formula for 
this model. 
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In order to avoid possible confusion, we call the maximally skewed Levy process 
with a = 1/2, j3 = ±1, the Levy-Gauss Process (LGP). By doing so we emphasize 
the so-called duality relation between the Levy process with a = 1/2, /3 = ±1, and 
the standard Gaussian process with a = 2, see, e.g., [122]. 

Maximally skewed stable processes possess are scaling invariant, which allows 
one to write their Green's functions in a simple form. Interestingly, a maximally 
skewed tempered stable process (c + = and k_ > 0) can be transformed into a 
maximally skewed stable process (c+ = and k_ =0). This result, which is listed 
in Appendix B, allows us to represent the Green's function for maximally skewed 
tempered stable process in closed form. It also allows us to derive closed-form call 
option prices for certain TSPs, which shall be useful later on for testing purposes. 
For instance, we have: 

Proposition 4.3. Consider the price of a T -maturity, K-strike call option in an 
exponential tempered Levy-Gauss process with a = 1/2 and c + = 0. Set l = \T23t , 
i\ = ^/kZl, t2 = yjn- + It, and k = ln(K/F(t)). Further, assume that k < 12 — i\, 
and set 

2l 2 

v = . 

O2 - i-i-k) 

We then have 

(4-14) ^p=C GL (r,k), 

C GL (t, k) = e 12 D (v, 2i 2 ) - e k+Ll D (v, 2u x ) , 

where D (v, k) is the symmetrized version of the Black-Scholes price given by 

(4.15) D(«,fc) = e-* (l-C BS {v,k)) = e -**(-d+) + e*$(d_). 

Proof. See Appendix B. □ 

Similar formulas, albeit in a somewhat different context are given in [54]. 

Remark 4.2. When a = 1/2, c + — 0, the price process F(T) is always bounded 
from above by 

F max {t,T)=F(ty*- l \ 
A similar bound exists for all a < 1. Paradoxically, when a > 1 and c + = 0, F(T) 
can take any value in (0, 00), despite all jumps being strictly negative. 

4.2. Normal Inverse Gaussian Processes. As we have seen earlier, in some 
cases Levy processes (such as a tempered Levy-Gauss process) are localized on a 
semi-axis. Processes of that nature can be used in order to create new interesting 
processes by changing time for the regular Brownian motion. Below we discuss some 
examples along these lines. A useful aspect of this approach is that it generates 
Levy processes with PDFs known in closed form. One of the well-known examples 
is the so-called Variance Gamma process (VGP), which is a special case of a TSP 
with a = 0, see, e.g., [82]. Another popular choice is the so-called Normal Inverse 
Gaussian process (NIGP), see, e.g., [12], [42], which we consider in some detail. 

The Inverse Gaussian (IG) process describes the density of the hitting time r 
of a level x\t by a standard Brownian motion with volatility o\ and drift x\. Our 
choice of parameters ensures that E (r) = t. It is easy to verify that IGP as a 
TSP with a = 1/2, /? = 1, k + = k\/2<j\. It is clear that IGPs and TLGPs arc 
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closely related and can be transformed into each other via shift and reflection. IG 
process can be used as a subordinator in order to build the N1GP out of a standard 
Brownian motion. 

The NIGP can be obtained by time changing of a standard Brownian motion 
with volatility a and drift — er 2 /2. The drift is chosen in a way ensuring that the 
corresponding NIG is a martingale. The corresponding subordinator is distributed 
as a hitting time of a level x\t by an independent BM with volatility u\ and drift 
xi . The convolution of these two processes yields 

(4.16) f NIG (t,x) = **ff~trT gl (W* 2 + *W) , 

where K\ (.) is the modified Bessel function, and Co = yjx 2 + 1/4, k— K\ja\G. 
The corresponding Levy exponent and density have the form 



^ NIG (u) = <T 2 X (x - ^X 2 + U (u + l)} 



V NIG (dx) = ( lim 1^1) dx = 1 K x (Co \x\) dx, x ? 0. 

Although NIGPs are not a special case of TSPs, they are closely related to TSPs 
with a = 1/2 and a = 1, due to the fact that 



Kl(Co\x\) ~ -J-f, K!{u\x\) ~ ( 1 + TT^TT ) 

v 1 l7 |xKow|a;| v 1 l; |x|-coo V 2w b V 8Cd\x\ > 



so that 



^ G(dx) „ ^/c^) „ e -(|-w>N). 



kl->o 7T \x\ kKoo ^/27r|a; 

In the spirit of formula (1.6), the price of a call option with log-strike k can be 
written in the form 

(4.17) ^^ = C»-( , l;i ), 
where, as usual, u = a 2 r, and 

(4.18) (v, *; *) = / (e /- e 7 Xi (* v^TW) 



Clearly, the normalized call price for a NIGP C NIG depends on one non-dimensional 
parameter. (Recall that C BS does not depend on any parameters.) In Section 5 we 
derive an alternative integral representation for C NIG , given by the Lewis-Lipton 
formula. We use one or the other of these two expressions as convenient. 

As before, it is easy to compute the annualized standard deviation of a NIGP. 
The corresponding expression has the form 

(4.19) stdev(X(l)) = —. 

K 

Remark 4.3. We note in passing that the so-called Generalized Inverse Gaussian 
( GIG ) processes can be viewed as a natural generalizations of IGPs. In general, the 
density of such a process is not known analytically, while its Levy density is given 
by a fairly complicated expression, so that GIG process is not easy (or necessary) 
to use in practice. It is not in the TSP class in any case. 
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4.3. Merton Processes. One popular process that is not related to the tempered 
stable class is the (pure-jump) Merton process (MP), with 

v M {dx) = -—e~ J ^dx. 

This process has finite activity, so we may use (3.15) to establish ip(u). We get 

[e iux - 1] u(dx) = A - l) , 

so, for a pure-jump MP, 

(4.20) i){u) = A (e^-W^ 2 _ fi _ e ^+^ 2 ) iu = \ (jw-W*' _ {^ +liu ^ 

where 7 = A (1 — e q ) , q = + ^i] 2 . The corresponding PDF can be calculated 
explicitly, 

n=0 



I 6 (x — "ft) , n = 0, 

fn(t,x) = { _ (»-7t-n„) 2 

1 1 -e , n > 0, 



/2-Knri 



It can be shown easily that 

v M (dx) = flim IS^T] dx = \f x (0,x) = — ^e -1 ^, a: ^ 0. 

A simple calculation yields the following expression from [94] for the price of a call 
option 

(4-21) ^ =C {v,k;\,v,r]), 

OO 

C M («, fc; /x, 77) = ^ c„C ss (nr/ 2 , Z - ng) , 

n=0 

where 

C = e- e9t, (e^)"/n!, i = fc-(l-e«)u, v = At. 

Analysis of MPs with diffusion component is very similar to the one performed 
above. The corresponding characteristic function has the form 

i){u) = -\o 2 u (u + i) + \ [e^-W^ 2 - l) + liu, 

while the price of a call option can be written as follows 

C M (t;T,ir; A, At, 77) _ M „ 
^ =C {v,k;a,fi,r]) , 

00 

C M (w, fc; <7, /i, 7?) = ^ c„C BS (<t 2 u + nrj 2 ,l - nq) , 

n=0 

where c„, I have the same form as before, and a = a j\f\. 
The annualized standard deviation of a MP has the form 

(4.22) stdev (X (1)) = (a 2 + A ( f i 2 + V 2 )) 1/2 . 
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Once again, this formula allows us to get a rough estimate of the magnitude of the 
implied volatility of ATM options on MPs. 

5. The Lewis-Lipton Option Price Formula and its Implications 

5.1. Exponential Levy Processes. The key formula allowing one to analyze 
option prices for Levy processes is known as the Lewis-Lipton (LL) formula; it has 
been independently proposed by Lewis [72] and Lipton [76] . This formula is based 
on the Fourier transform of an appropriately modified payoff of the call option. A 
complementary method is proposed in [27]; additional information can be found in 
[68]. 

Proposition 5.1. The normalized price of a call option written on an underlying 
driven by an exponential Levy process with Levy-Khinchine exponent tp has the form 

(5.1) C (r,fc) = l-^- /" ^4e- k (^)du. 
Here 

(5.2) E (t, u) = cxp (to (it)) = cxp ^v (u) - ( u ) 

(5.3) v(u) = (u- , v (u) = 4> (u - £j , Q{u)=u 2 + ^. 
andip is given by (3.4-). 

Proof. Sec [71], [72], [75], [76]. □ 
As a corollary we have: 

Corollary 5.2. The derivatives with respect to k of the normalized price of a call 
option are 

,5.5, c l(T ,t)-4£^.-K-.)(- ta+ i)V 

Applying this result to the processes introduced in Section 4 we easily get the 
following lemma. 

Lemma 5.3. For TSPs, NIGPs, and MPs, the corresponding E(t,u) have the 
form 

E ts {t,u) =exp ^-l(T 2 Q(u) + J2 a s ( /s «~ s ( iu +^)) +7 ( iu+ ^) +6 
E NIG (t, u) = exp (t<t 2 k (x - ^fuo 2 + it 2 ) ) , 
E M (r, u) = exp (r (-\<? 2 Q («) + A (e< iu+ ^ Q ^ - 1 + (1 - e") (iu + 
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Remark 5.1. It is worth mentioning that E TS (r, u) and E NIG (r, u) rapidly decay 
when \u\ — > oo, so that computation of the integrals (5.1), (5.4), and (5.5) is 
straightforward. When a ^ 0, E M (r, u) is rapidly decaying as well. However, 
when a = 0, the situation is more difficult. The corresponding integrals (5.1), (5.4) 
for this case can be computed efficiently by using 

E M ^ (t, u) _ cxp (At (-! + (!- 

qJu) oJu) ' 

as a control variate. 

Remark 5.2. In the case of a BSM diffusion process, (5.1) yields 



1 



^ vQ(u) — k(iu— |) 



C BS (v, k) = 1 - -f- / du, 

1 ' InJ-oo Q(«) 

so that 

5.6) — / — = $ — - -y/v )+e k <£>[ 1 =--< s /v\. 

In i/ie limiting case v — 0, we get 

e -fe(»«-|) 



1 A 



d?i = e fe . 



Q(u) 

These useful formulas shall be used repeatedly in what follows. 

The usefulness of the LL formula becomes apparent when one wants to study 
the asymptotic behavior of the call price and its derivatives (the Greeks) in the 
limiting cases of r — > 0, t — > oo, \k\ — > oo, since it allows one to use enormous body 
of work dedicated to the asymptotics of integrals depending on large and small 
parameters. This is done in the remainder of the paper where it is shown how to 
apply the saddlepoint method, the high frequency Fourier integrals estimates, and 
other tricks to the problem at hand. We return to these asymptotics in Sections 7, 
8, 9. While this is not the focus of this paper, the LL formula can also be used to 
study the small jump asymptotics, as is briefly shown in the next section. 

5.2. Small jumps asymptotics. When the jump component of a Levy process is 
small compared to its diffusion component, by expanding E (t,u) in (5.2), the call 
price can be written in the form 

1 r°° p — jC 2 t<2(m) — k(iu— 5) 

(5.8) C(r,fc) = l-— / — du 

' y ' 2W-oo Q(«) 

1 r°° e -5 <T2r Q(")- fe (" i -5) 

-T— pr-T-T V (u)du+..., 

2tt 7_oo Q (u) 

where vq (u) is given in (5.3), provided that the second integral converges. This 
expression is particularly useful for the qualitative study of perturbations of the flat 
volatility surface caused by jumps. Indeed, by comparing (5.8) with the expansion 
of the BS formula around Oi mp — a and matching terms, the implied volatility 
surface can be represented in the form 

o-imp (r, k) = a + cti (t, fc) + .... 
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where o\ is of the same order of magnitude as vq, and is given by the following 
expression 

a x (r, k) = t 1 ' 2 — '=Lr / — — v (u) du 

V27T J-oo Q( u ) 



-§(- 1/2 «+^ 



-t>o (u) rfu. 



Q(u) 

A similar formula, albeit derived in a much more complicated way, is given in [90]. 

5.3. Quadratic Volatility Processes. On rare occasions, formulas similar to the 
LL formula can be used for spatially inhomogeneous processes as well. The best- 
known example are the so-called quadratic volatility processes (QVPs), which we 
shall briefly describe below. The reader is referred to [76] and [3] for further details. 
The reason why these processes are considered alongside ELPs is due to the fact 
that after appropriate transforms they can be made translationally invariant, see, 
e.g., [75], [76], [26]. 

When the volatility is quadratic, (including the limiting case when it is linear), 
the dynamics of the corresponding underlier is driven by the following local volatility 

sde' 

dF(t) = a? oc (F(t))dW(t), F(p) = F , 

where 

<jf oc (F) = aF 2 + bF + c, a > 0. 
Consider the following quadratic equation 
(5.9) <(F)=0. 

If we want to ensure that af^ c (F) > for < F < oo, we have to consider two 
possibilities: (A) (5.9) negative real has roots; (B) (5.9) has complex roots. For 
brevity, we concentrate on the second possibility, which is most often what real 
market data dictates, and write 

aL(^) = a(V-p) 2 + q 2 ), q >0. 
It turns out that in the case in question the following proposition holds. 

Proposition 5.4. In the case when volatility is quadratic with complex roots, we 
can represent the price of a call option using an eigenfunction expansion represen- 
tation in the form 

(5.10) 

C(t;r,K) t 1 



F t (p sin (A Ft ) + qcos (A F J) 

R(h) 



~ e -lvR(ki) 

R(h) (C" sin ( 2k i X K) - 2fc ; C c cos (2k[X K j) sin (2k t X F ) , 
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where F t = F (t), v = 4a 2 q 2 r, k t = Trl/2X ca , and 
X Ft = ^arctan ^ — — + arctan 



X K = I 



, K-p\ fp 
arctan + arctan 



1 J \q 



, 7T fp 

-X-oo = t: + arctan 



2 \q 
C = t^f- ((P 2 + q 2 - Kp) cos + Kg sin (X K )) , 

C = ((P 2 + q 2 - tfp) sin - Kqcos (X K )) . 

froo/. Sec [76], [3]. □ 

Remark 5.3. Expression (5.10) is more compact, but equivalent, to the one given 
in [3]. As in [3], solution based on the method of images is also possible. As usual, 
the eigenfunction expansion based solution and the method of images based solution 
should be used when r — > oo and r — > 0, respectively. 

Remark 5.4. It is worth noting that drift- free processes with quadratic volatility 
are not martingales, but rather supermatingales, see., e.g., [3]. ,4s a consequence, 
the solution of the corresponding pricing problem is not unique and a proper one has 
to be chosen. Such a choice is implicitly done above. Non-uniqueness also means 
that there are (many) non-zero solutions of the pricing problem with zero initial 
and boundary conditions. To demonstrate, assume for brevity that p = q = 0, so 
that 

a^ oc (F) = aF 2 , a>0. 
The corresponding homogeneous pricing problem can be represented as follows: 

V T - \a 2 F A V FF = 0, 

V(0,F) = 0, V(t,0)=0. 

It can be shown that the generic non-trivial solution of the above problem has the 
form 

1 r e -{^F 2 {r-r Or 

V (t, F-h) = - / h (r') dr'. 



a Jn >2n (r - r'f 



In particular, when h (r) = 1, we have 

V(t,F;1) = 2F$ 



aF^r 



5.4. Heston Stochastic Volatility Processes. Another important case when 
the LL formula can be used successfully is the so-called Heston model 9 . Heston 



We note in passing that many of our results arc applicable for other cases, for instance, for 
the Stein-Stein stochastic volatility processes. 
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stochastic volatility processes (HSVPs) are governed by a system of SDEs of the 
form 

dF (t) = \fwj€)F (t) dW^ (t) , F{0)=F , 

dzu (t)=K(6-w (t)) dt + Eyjvo (t)dW (ru) (t) , zu(0) = zu 0l 

where dW^ (t) (t) = pdt, see [59]. According to [75], [76] for the Heston 

model we can represent the price of a call option as follows. 

Proposition 5.5. The normalized price of a call written on an underlying driven 
by a square-root stochastic volatility process has the form 

1 f°° E(t,u) k ( iu _i) 



(5.11) 






where 






(5.12) 


E(t, 


u) 




A(t, 


u) 




B(r, 


u) 




C(r, 


u) 



C (t, k) = 1 - — / v .'/ e -H'"-iJdu, 



e A(T,u)-B(r,u)vj Q(u) _ e C(r,u) 



K9 ( F+ («) r + 2 In ( T ~ ^ + T+ ^ ° XP iy ~ Z ^ ^ 



2Z («) 
1 — exp (— Z (u) t) 



T- (u) + F + (u) exp (— Z (u) r) ' 



and 

(5.13) J-±(«) = ± (pel iu+l ) -k) +Z(u), 



or, equivalently, 



Z{u) = \l (ps[iu+l) -k) +e 2 Q{u), 



F± ( u ) = ± (peiu — k) + Z (u) , 



Z (u) = \J {peiu - kf + e 2 Q (u), 

where k = k — pe/2. 

Remark 5.5. The martingale condition, which is easy to verify, reads 

Remark 5.6. We note in passing that C (r, u) can be represented in the form 

C (t, u) = OCi (t, u) + tu C 2 (r, u) , 

which emphasizes the contributions of average and instantaneous variance, respec- 
tively. 

It is clear that C (t,u) in Proposition 5.5 is not a linear function of r, so that 
the Heston process is not a Levy process. However, in the limits of r — >■ oo and 
r — > it can be viewed as such. 
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Proposition 5.6. For t — > oo we can represent C (t,u) as follows 



(5.14) C(t,u) = - 
For t — )■ we have 

(5.15) c(t, V -)=- 



k9T+ (u) t 2k9 



iw Q v sinh ( e ^ i ) 



In 1 



2Z (it) 



G7 



7"+ («) + O 



+ 



re sinh 



v 2 



2 



sinh p- 
pe 2 sinh ( 
arctan (p/p). 




+ 0(t), 



where p = \fl — p 2 , and (j> - 

Proof. Straightforward but tedious calculation leads to (5.14). To derive (5.15), we 
use for inspiration the well-known duality between the Brownian motions W (t) 
and tW (1/t), introduce a new variable v = tu and then the formula follows. □ 

It is easy to see from formula (5.14) that a Heston process with zero correlation 
is asymptotically equivalent to an NIGP. Naturally, in the long-time limit, to the 
leading order C does not depend on ru n . Thus, to the leading order a Heston process 
can be viewed as a Levy process but with time inverted. Also notice that in the 
short-time limit C does not depend on 9 to the leading order. 

A useful survey of modern approaches to option pricing in the Heston framework 
is given in [121]. 

5.5. Calibration to Market Data. For later tests we shall need concrete parametriza- 
tions of the various processes introduced in this paper. For this purpose, let us 
attempt to calibrate to the set of market data given in Table 1. Since we arc only 
interested in the case when parameters are constant in time, it is not possible to 
calibrate any of the processes of interest to the entire set of market quotes. Instead, 
we shall choose a representative maturity, r = 2y, say, and calibrate parameters 
to the set of five option volatilities. This will allow us to identify proper order of 
magnitude for these parameters. We consider six representative ELPs and present 
the corresponding calibrated parameters in Table 3. 



TSP1 


a = 0.66 


a = 


0.07% 


c+ 


= 0.1305 




0.0615 


K+ 


= 6.5022 


k_ = 3.0888 


TSP2 


a = 1.50 


a = 


0.01% 


c+ 


= 0.0069 




0.0063 


K+ 


= 1.9320 


k_ = 0.4087 


MP 


A = 35.33% 


1* = 


-0.0318 


V 


= 0.2023 












HSVP 


k = 2.2707 


9 = 


0.0225 


e 


= 0.6200 


P=- 


-0.0541 




= 0.01374 




NIGP 


a = 14.90% 


X = 


= 3.20 
















QVP 


a = 1.322 


P = 


0.967 


q 


= 0.301 













Table 3. Parameters for six representative processes calibrated 
to market quotes at t = 2y. 



The quality of the calibration for two TSPs, MPs and HSVPs, which is fairly 
high, is shown in Figure 3(a). Although for NIGPs and QVPs do not have enough 
parameters to ensure successful calibration to the market, our choice of parameters 
generates satisfactory (but not perfect) fits shown in Figure 3(b). 
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(a) TSP,MP,HSVP 




(b) NIGP,QVP 



Figure 3. a (A) calibrated to five market quotes (see Table 1) for 
four representative processes. The corresponding maturity is 2y. 
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6. asymptotics of the option price in the black-scholes-merton 

Framework 

In Sections 7, 8, 9 we shall use the LL formula of Section 5 to establish a variety 
of asymptotic expressions for option prices and implied volatilities. Before this, 
however, we take a brief detour into the analysis of the limiting behavior of BS 
prices when variances v and strikes k are either large or small. This analysis is 
required later in order to convert option price asymptotics into implied volatility 
asymptotics. 

As was mentioned earlier, the undiscounted price of a call option can be written 
in the form C BS (v, k), where v is the annualized variance, and k is the log strike. 
The well-known relations 



(6.1) 



$ (x) = 1 - 4> (-x) , $ (x) 



X 



1 



1 



allow us to simplify the above formula in various asymptotic regimes. Specifically, 
we are interested in the following three cases: A) v — > oo, k = kv, k is fixed 
(long-time asymptotics); (B) v — > 0, k is fixed (short-time asymptotics); (C) v = 
v\k\, J-oo,0<t)<2 (wing asymptotics). Specifically we have the following 
proposition: 

Proposition 6.1. In case (A) we have 
(6.2) 



•BS 



(l -e*") !{(<_■} + 1 {-;<i<;} 



v R(k) y R(k) 
M } + (l-^/ 2 )l {fc= _ J} 

fi-l/2)t</2 



k^±\ 



(1-^)+-. 



2 " 



k — ± 2 , 



R (fc) = -Q (ifc) = k 2 



where 
(6.3) 

In case (B) we have 
(6.4) C BS (v, k) = (1 - 

Finally, in case (C) we have 
(6.5) 

Q BS {v\k\,k) = {l-e k ) + + - VV 



1 

4' 



6 2 w 3/2 
k 2 



\k\ 



1 



R{l/v)^i\k\ 



1 - 



3/v 



r,2 



{R{l/v)Yv\k\, 



+ .... 



Proof. Simple but slightly tedious application of (6.1) to BS formula (1.6), see, e.g., 
[117], [47], [52]. □ 

Remark 6.1. For case (A), it is clear that cases k = ±1/2 require separate treat- 
ment, with k = 1/2 being particularly important since in order to get A ~ 1/2, we 
have to have k — 1/2 ~ \ j^fv, in which case 



A ~ $ \k- 



1 
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Remark 6.2. For future reference, we wish to generalize formulas (6.2) and (6.4) 
and analyze the asymptotics ofC BS (v) v, kvj , with g(v) = (l + a.\jv + a 2 /v 2 ), 

k is fixed and v — > oo, and C BS (f (v) v, k) with f (v) = (l + a\V + a 2 v 2 ) , k ^ is 
fixed and v — > 0. Straightforward computation yields 

(6.6) 

(g («) v, kv) = l {k< _ h} (l - e-) + l { _ h<k<h} + e^«« " j_ ' 




I 1 -II — txn — -tt I |3+^rv I +^ + ^(«?-« 2 ) 



and 



(6.7) C BS (f(«) w ,fc) = (l- e fc ) + + ^A^ 



fc 2 



1- * (1-12^+4 (a 2 - a 2 )fc 2 )v 

These formulas are used below for studying the asymptotic behavior of the implied 
volatility in the long-time and short-time limits, respectively. A more complicated 
expression, which is equivalent to (6.6) is given by [47]. 

7. Long-time Asymptotics 

7.1. General Remarks. For the purpose of establishing long-time implied volatil- 
ity asymptotics via the LL formula, let us briefly introduce the so-called saddlepoint 
method. 

Saddlepoint approximation is a method for computing integrals of the form 
(7.1) g(r) = ^J f(z)e^dz 

when t — > +oo. Here £ is a contour in the complex plane, and the amplitude and 
phase functions f (z) ,S (z) are analytic is a domain V containing £. The extremal 
points of S, i.e. zeroes of S' , are called the saddlepoints of S. Under reasonable 
conditions, the contribution to g (t) from a vicinity of a saddlepoint zq, where 
S" (z ) ^ 0, is given by 

<7 - 2) s >--° iT) = 7^m erSM,M{1+oiT ~' )h 

see, e.g., [43], [80], [63], It is clear that the main contribution to the integral comes 
from the saddlepoints where Re [S] attains its absolute maximum. The second-order 
approximation has the form 
(7.3) 



, , 1 j (AS)' , 5 (*"'(*>)) Qgo) ■ . , , 



-2\ 
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Saddlepoint approximation has been successfully used by many authors for a 
variety of financial applications, see, e.g., [104], [71], [28], [86], among others. 



7.2. BS Asymptotics via the Saddlepoint Approximation. To motivate sub- 
sequent developments, let us briefly derive some large v asymptotic results for the 
BS case. We are specifically interested in the case when v — > oo, k = kv, k is fixed. 
We have 

C BS v, kv) = 1 - ±- / nl ) du = 1 - ^ / e^ k ) f (u) du, 

V ' 27T J_ x Q (U) 27T J V ' 



where 



(u,Jej = ~^Q(u) - k (iu- , f(u) 



Q( u y 

The integrand is meromorphic in the entire complex plain where it has two simple 
poles located at the points 

u± = ± 1 -, s(^,kj=k, s(- l -,kj=o. 

We calculate the location of the saddlepoint by solving the equation 

S' (u, kj = -u -ik = 0, 

so that 

«. = -»&, S" («.) = -!, S(u t ) = -Uk-^j , /(«.) = 



i? 



(t) 



and i?(fc) is given by (6.3). In order to apply the saddlepoint approximation, we 
need to transform the contour of integration in such a way that it passes through 
the saddlepoint, so that its immediate vicinity dominates the entire integral. We 
can achieve this goal by parallel shift of the contour of integration from the real 
axis to the contour given by 

Im (u) = —k. 

In the process of doing so, three possibilities might occur: (A) the lower pole tt_ 
is crossed (OTM case); (B) no poles are crossed (near- ATM case); (C) the upper 
pole U- is crossed (ITM case). When poles are crossed, their contributions have to 
be computed via the Cauchy formula: 

n+ = e kv , n_ = 1. 

The contribution of the saddlepoint has the form (7.3). Summing up all these 
contributions, we obtain formula (6.2) of Section 6, where it is derived by elementary 
means. 



7.3. Exponential Levy Processes. We can now consider more general LPs. 
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7.3.1. Generic Exponential Levy Processes. For a generic ELP the LL formula 
yields 



where 



1 r 

C(r,fc) = l-— / e^f(u,k)du, 

^ J -co 



1 p~ k { 1 
S(u) = v (u)--a 2 Q(u), f(u,k) ' 



2 " w ' J v ' ' Q(u) ' 
where (w) is given by (5.3). Thus, for large r the integral of interest has the form 
(7.1), and can be analyzed via formulas (7.2) or (7.3). However, it is more natural 
to assume that k ~ kr, so that the strike moves deeper in or out of the money 
when maturity increases (this case is necessary to consider in order to study the 
asymptotic behavior of volatility as a function of delta) . Under this assumption we 
have 



where 



C(r,fc) = l-i-| ° e TS (^ k )f(u)d Ul 
S (u, k) = v (u) - t,o- 2 Q (u) -k^iu-^J , / (u) 



Q{uY 

Consider now the general case of ELPs. Since we are dealing with shifted charac- 
teristic functions, we can use symmetry and prove that the saddlepoint of interest is 
located on the imaginary axis. Accordingly, on the interval of analyticity of S (u), 
given by the inequalities Y_ < Im (u) < Y + , we can define a real- valued function 
Soi (y, k) of real arguments y, k (since on the imaginary axis the value of S (u, k) 
is real, see, e.g., [71], [81]) as follows 



Soi (y, k) = s (y) + k (y+\^j 



where 

(7.4) Z a (y)=v (iy) + ^a 2 R(y), 

and R (y) is given by (6.3). It is clear that we can find the location of the saddlepoint 
by solving the following equation 

(7.5) S (y) + k = 0. 

Since, in general, this equation cannot be solved analytically, we prefer a parametric 
approach, by expressing k in terms of y, rather than the other way around. This 
approach leads to the following proposition. 

Proposition 7.1. Consider y e (— Y + , Y_) and define 5 (y) , Si (y) andE i^± (y) 
as follows 

(7.6) So (y) = vo (iy) + l -a 2 R (y) , E 1 (y) = S (y) , 

1 



S i,± (y) = S {y) - Si (y) [y± 2 

Then for k = —Si (y) the corresponding a imp (r, rfc) can be written in the form 
(7.7) a lmp (r, rk) = (a (,) + ^ + ^ + ...) ^ , 
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so that a i m p (r, rfc) can be represented as follows 

(7.8) a imp (r,rk)=bo(y)+ b ^ + b ^ + .... 
where 

(7.9) a (y) = (^sign (y + yJ-2~oi,+ (y) - sign (y - ^/-2S 01 ,_ (y) 

(7.10) o (») = <4 /2 (y) , 

and higher order coefficients Oj (y) and 6j (y) /iawe ifte /orm (D.6), (D.7). 

Proof. See Appendix D.l. □ 

Fundamentally, Proposition 7.1 is derived by comparing the asymptotic expan- 
sion obtained via the saddlepoint method with the asymptotic expansion for BS 
price given by (6.6), and balancing terms. Below we use the following notation 

(7.H) o- imp , i (y) = J2 b ^, 

3=0 T 

for partial sums of an asymptotic series (7.8), and similarly in other instances. 

Remark 7.1. While in our approach Proposition 7.1 follows from asymptotic for- 
mulas for complex integrals, it can be obtained via the large deviation principle as 
well, see, e.g., [118], [119], [37]. 

7.3.2. Specific Exponential Levy Processes. For specific ELPs, such as TSPs, NIGPs, 
MPs, etc., the corresponding formulas can be made explicit. 

Proposition 7.2. For TSPs, NIGPs, and MPs Proposition 7.1 holds provided that 
the corresponding S (y) are defined as follows: 

1 . ( ( \\\ a ( 1 



2 ■ ^ - v - v" ■>)) '{" 2 



(7.12) ~l s {y)= L -a*R{y) + Y, 

S— i 

(y) = a 2 x (x - ^Q 2 - y 2 ) 



~NIG 
"0 



S M (») = \o 2 R (y) + A (e-«(v-h)+4"(v) - 1 - (1 - e *) ^ - | 

for TSPs, y e (-K+ + |,K_ + i); /or ATCPs, y e (-w, u>); /or MPs, y e (-oo, oo) . 

Proof. Proof is straightforward. For further details, see Appendices D.2, D.3, D.4, 
respectively. □ 

The quality of the saddlepoint approximation in the limit of infinite maturity 
for representative TSP, NIGP, and MP is illustrated in Figures 4, 5, 6, and 7, 
respectively. The relevant parameters are calibrated to the market. It is clear that 
for the calibrated TSP with a = 2/3, NIGP, and MP saddlepoint approximation is 
very good, while for the calibrated TSP with a = 3/2 it fails. 
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Figure 4. Comparison of the exact and asymptotic expressions 
for <7i mp for the calibrated TSP for different maturities, a = 2/3. 
Here and in Figures 5, 6, and 7, "exact" denotes the implied volatil- 
ity calculated by virtue of the LL formula, while a Si, a A, aJl are 
given by (7.11). 
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Figure 5. Comparison of the exact and asymptotic expressions 
for Gimp for the calibrated TSP for different maturities, a = 3/2. 
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Figure 6. Comparison of the exact and asymptotic expressions 
for Gimr, for the calibrated NIGP for different maturities. 
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Figure 7. Comparison of the exact and asymptotic expressions 
for aimr, for the calibrated MP for different maturities. 



EXPONENTIAL LEVY PROCESSES 



37 



Remark 7.2. For some special cases the above formulae can be made more specific. 
For the general maximally skewed TSP we can solve (7.5) explicitly and avoid using 
the parametric representation. Specifically, we have 



s {y) = a 



Si (y) - 

where (, r\ are given by (4-5), so that 




1 


-a- 




2 




r- 

V a 




(k) - sign (k+ - k) yJ-2Eo!,- (k). 



These formulas show that for TLGP the ATM implied volatility approaches from 
below the following level 



^(0) = 23/V/ 4 (v^TT-^ 3/2 



We compare asymptotic and exact implied volatilities for TLGP in Figures 8 
(a), (b) below. These Figures show that, similarly to the general case of TSPs with 
< a < 1, for TLGPs the saddlepoint approximation is acceptable. 
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exact a_0 — o_l o_2 







(b) T=15y 

Figure 8. Comparison of the exact and asymptotic expressions 
for crimp for the GLP with d = 0.0075, k = 1.0 for different ma- 
turities. Here <Ji m p,exact denotes the exact volatility calculated by 
virtue of (4.14), while <J imPii are given by (7.11). 
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Remark 7.3. For NIGP we can easily invert the equation S (y) + k = 0, and 
avoid using a parametric representation. We have 

yo{k)= - F =dL=== , S i,± (k) = a 2 x 2 - usj ct 4 x 2 + k 2 ± \k, k± — T~7r, 

V cr 4 x 2 + k 2 z z 

^ G co (k) = sign (k_ - k) ^/-2S 01i+ (k) - sign (k + - k) ^/-2S 01 ,_ (fc). 

Remark 7.4. 5y using (4-18) which provides the call price for a NIGP explicitly, 
we can calculate long-time asymptotics for call prices directly, see Appendix D. 3 for 
details. The asymptotics obtained via two complementary methods naturally agree. 

7.4. Heston Stochastic Volatility Processes. Equation (5.12) shows that for 
the Heston model the LL exponent is not proportional to time. However, it is 
proportional to time to the leading order when r — > oo. Formal expansion in 
powers of exp (— tZ (uj) yields 

E(T,U,k) _ e Mr,u)-B(r,u)^ Q(u)-rk{iu-l) 



Q (u) Q (u) 



r.e T ^>f(u), 



where 



(7.13) S(u,k)=-^T + (u)-k(^iu-^j 



/(«) 



2k6/s 2 ' 



and J+ (u) , Z (u) are given by (5.13). As before, we can use saddlepoint method to 
obtain the asymptotic of the LL integral. It is easy to check that the corresponding 
saddlepoint has to be purely imaginary, so that we can proceed as before. 

Proposition 7.3. Consider y e (— Y + ,Y_), where 

Upk + Jk 2 + \p 2 ) 

(7.14) Y± = ^ V j 2 

and k = k/e, and define 

k8 

(7.15) So(») = — (py + k-<;(y)), 

^ , x K0 ( -p 2 y + pk 



s i,± (y) = s (y) - Si (y) (y±\ 
where 

(7.16) ? (y) = ^-p 2 y 2 + 2pky + k 2 + 1/4. 

Then for k = —Si (y) ifte corresponding cr imp (r,rfc) can 6e written in the form 

a mp (r, rfc) = (a (y) + ^ + ^ = b (y) + ^ + 

where the leading order coefficients a$ (y) and bo(y) have the form (7.9), (7.10), 
and higher order coefficients ai (y) and bi (y) have the form (D.6), (D.7). 
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Proof. The proof is similar to the one of Proposition 7.1 and is omitted for brevity. 

□ 

As expected, the leading-order term in the volatility expansion does not depend 
on wq. The quality of the above approximation is illustrated in Figures 9 (a), (b). 
These Figures show that for HSVPs the saddlepoint approximation works well. 




(a) T=2y 




(b) T=10y 



Figure 9. Comparison of the exact and asymptotic expressions 
for (Ti m p for the calibrated HSVP for different maturities. Here 
"exact" denotes the implied volatility calculated by virtue of the 
LL formula, while cr_0, cr_l (7.11). 
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Remark 7.5. For HSVPs we can easily invert the equation Sq (y) + k = 0, and 
avoid using the parametric representation. A simple calculation performed in Ap- 
pendix D.5 yields 

(7.17) 

o-i m p,o (&) = sign (y* (k) + ^J-2E 0h+ (fc) - sign (y* (k) - tJ-2Z 01 _ (k), 
? n (0 y^ 2 - p 2 (« 2 - p2) / {p 2 + F) 



pn - sigr 

y* (k) = — 



P 2 



kO 1 - 

Hoi,± (fc) = Tl (if + pi) y* (k) + (l-p)k)± -k, 

where I = p + ek/nO. 

Additional information is given in [51], [62]. 

8. Short-time Asymptotics 

8.1. General Remarks. From the results in Proposition 5.1, and Corollary 5.2, 
we see that in order to analyze the asymptotic behavior of the ATM option price 
(k = 0) and its derivatives for r — > 0, we need to study the following integrals 

(8.1) l(r) = ± r {l -^ T : U)) d u = 1 Re ( f ^j^du) , 

V ' 2ttJ_ 00 Q(u) tt \J Q(u) J ' 

- (r) = ^ f = - 1 Im( ^ ^udu 

V ' 2ttJ_ 00 Q(u) tt Uo Q(u) 

If 00 If/" 00 1 

n (r) = — / E (t, u) du = - Rc < / £ (r, u) du > , 

27r J-oo t Uo J 

where 

/ <7 2 T 

E (t, u) = exp I ru (u) —Q (u) 

The relation between C (t, 0) , Cfe (t, 0) , C^^ (r, 0) and I (r) , m (t) , n (r) is straight- 
forward. 

Proposition 8.1. C (t, 0) , (r, 0) , (r, 0) can &e expressed in terms of I (r) , m (r) 
as follows 

(8.2) C(r,0) = I(r), 

C fe (r,0)=m(r) + i(I(r)-l), 
Q fe (r,0)=n(r)+m(r) + i([(T)-l). 
Proof. Straightforward computation. □ 
We assume that the implied volatility can be expanded in powers of k: 

, fen 

(8-3) a lmp (T,k) = 2^x„(t) 

n>0 

and demonstrate how to express Xo: Xi> Xi i n terms of [, m, n. 



42 



LEIF ANDERSEN AND ALEXANDER LIPTON 



Proposition 8.2. In the LL framework we get 
Xo(r)~i(r) ( 1+1 ^r) -~ 1/2 > 
Xi(r)~A(r) (*+^) ^ 

X 2 (r) ~ (- -L + ?M + J | (T) (T) + r" V, 

w/iere 

(8.4) I (r) = \/27r[ (r) , m (r) = v^m (r) , ft (r) = v^n (r) . 

Proof. See Appendix C. □ 

8.2. Exponential Levy Processes. 

8.2.1. Non-ATM Options on Tempered Stable Processes, Simple Heuristics. Using 
the LL formula, it is straightforward to develop small-time asymptotics for non- 
ATM options on TSPs. The argument goes as follows (see, e.g., [70] among others). 

Proposition 8.3. Assume that a = and consider 

i r°° e TV °( u ^~ k ( iu ~^ 

C(r,fe) = l-— / —- du, fc^O, 

where vq (u) is given by (5.3). If a € (0, 1), then 

d 1 f°° e~ k { iu ~^) 



dr 



1 f°° e ~ v ,u ~i) 
C (0, k) ~ -— / («) du, 

2?r J^co Q (u) 



w/iere the integral clearly converges. If a £ (1,2), i/ien integration by parts (which 
is possible because k ^ 0) shows that: 



CMJ-^.-K-HteMY*.. 



where the integral clearly converges. Accordingly, in both cases, the time value of 
the call option SC (t, fc) is linear in time, 

SC (t, fc) =c(fc)r + o(r). 

Proof. See Appendix E.l. □ 

Broadly speaking, this result means the time decay of the time value SC (r, fc) 
is much slower that in the BS case, which is of order exp (— fc 2 /2cr 2 r), see (1.6), 
(6.1). To compensate for that, the corresponding non-ATM implied volatility must 
explode with order 0(1/ y/r ln(l/r)), as was already mentioned in Section 1. 
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8.2.2. ATM Options on Tempered Stable Processes, Simple Heuristics . For ATM 
options we start with a simple heuristic argument. Consider a TSP with a = and 
a e (0, 1) and calculate the corresponding ATM asymptotics. From (3.10), we have 

(C(T,x + sy)-C(T,x))e-^y 



C T (t, x) = ~/C x (r, x) + ^c s 

Q=+ Jo 



s=± 

+ 



y 



l+a 



-dy, 



C(0,x) = (e*-1)- 
As we already know, the martingale condition yields 



s=± J 

where 



„ — 1_ jo y 



r°° (e s y - i)e- K °y , 

Ps = sc s / — dy > 0, s = ±. 

Jo V + 

In order to calculate C (r, 0) we discretize x on a grid x n = nh, —N < x < N, 
C n ( T ) = C(r, x n ), and use an explicit finite-difference scheme. Since we have an 
advection term, we have to differentiate two cases: (A) 7 > (p_ > p+), and (B) 
7 < (p- < p+). In case (A) our finite-difference scheme yields 

C (t) (Ci (0) - C (0)) (e h - 1) 

t h n 

In case (B) our scheme yields 

C (r) (Cp (0) - d (0)) 

= 7 r +P+ =P+- 

t h 

Note, in both cases differentiation is performed upstream to ensure stability of the 
scheme. We can combine the above formulae into one as follows 

(8.5) C (t, 0) = rmax(j) + ,p_) . 

We note that this formula coincides with (1.10) in Section 1. For finite variation 
TSPs, the ATM volatility therefore goes to zero (to ensure that C (t, 0) decays 
linearly), in marked contrast with non-ATM volatilities, which, as was noted in 
Section 8.2.1, explode for r — > 0. We emphasize that the technique above does not 
hold for a > 1, as the p± integrals diverge. We also note that the result gives us 
no information about the slope and convexity of the ATM implied volatility. To 
address these issues, we now process with a more detailed analysis. 

8.2.3. Main Result for Tempered Stable Processes. In view of the LL formula, in 
order to be able to evaluate the short-time asymptotic behavior of ai mp (r, 0) and 
its fc-derivatives dkOi mp (t, fc)| fe=0 , d\ui mv (r, fc)| fe _ Q , we need to evaluate the short- 
time limit of the three integrals in (8.1) with the shifted characteristic exponent of 
the form 

v (u) = ^ a s (k s - s (iu + ^ +7 (iu + + S. 

We have to distinguish four cases (A) a e (0, 1) , a = 0, (B) a e (1, 2) , a = 0, (C) 
a E (0, 1) , u > 0, (D) a G (1, 2) , a > 0. We calculate the corresponding integrals 
by scaling the independent variable u as appropriate. Specifically, we use: v = tu 
(case (A)); v = T x l a u (case (B)); v = t x I 2 u (cases (C) and (D)). This technique is 
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similar to the approach in [35], for example; in case (B), it has also been used by 
[116]. 

When considering the relevant integrals, we note that the dominated convergence 
theorem (DCT), which is the principal tool for studying the limiting behavior of 
integrals depending on parameters, cannot in all cases be applied to establish the 
limits for r — > directly, so great case must be taken to avoid faulty conclusions. 
To demonstrate this, let us just consider the simple example of establishing the 
limit of 

If l- e r T( m +2) 
2tt J u 2 + j 

— OO 

which is similar to l(r). Formula (5.7) shows that 

t(r) = l-e T7 ~ ~ -Tj-. 

Now, we wish to derive this formula asymptotically. We change variables v = tu, 
and represent I (r) in the form 

[(r)=r£(r), 

where 

„, , 1 f°° 1 - e? +i T" , 1 ri-e? cos hv) , 

We notice that at v = the corresponding integrand is unbounded when r — > 0, so 
that we cannot apply the DCT directly. However, we can split the above integral 
into two parts as follows 

£(r) = £g (r)+£r(r), 

where 

In order to evaluate £§ (r) we expand the integrand around w = and get, to the 
leading order in e, 

In order to evaluate £^° (r) we apply the DCT, which says that 

2?(t) -+ £f (0), 

T— >0 

since the corresponding integrands are uniformly bounded, yielding: 

7T.L V 2 £^0 2 



Thus, 



£(r)~=^±M = - 7 " 



which is the correct answer. 

Judicious usage of the technique above, leads to the following important result: 
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Proposition 8.4. Consider integrals I, m, n defined by formulas (8.1). Their as- 
ymptotic behavior depends on a and a. In cases (A)-(D) the corresponding asymp- 
totics has the form 







m 


n 


A 






9t T 


B 




°9Jt 


91 T 


C 




r (C) T l/2 
931 T 


1 | r (C)^(l-a)/2 


D 


^ 1/2 +C fW r (3-a)/2 
V2tt £ 


C rW T (2-a)/2 


1 + ciPM 1 - 01 )/ 2 

V2^ CTT i/2 ' 97 



where 









C<Xt 


A 


- (hr + S+g) 


- 5 sign (7) 


^ (7) 


B 


ir (1 - a') r Q ' cos (a'x) 




lT(l + a')r- a ' cos(a'x) 


C 


-(§7 + <5 + f?) 




I 2 (a-l)/2 r («±l^ tJ -(a+l, 


D 


_I 2 (a-3)/2 r (a f l^ (7 -(a-l) 


-i2( a - 2 )/ 2 r(f) 9 a- Q 


I 2 («-l)/2 r ( a±l) p(T -(a+l) 



.ffere a' = 1/a, a±,j,5, g are given by (4-5), 5jj (.) is the Dirac delta function, and 
p = (a+ + a_) cos < 0, g = — (a+ - a_) sin , r = \J p 2 + q 2 , x = arctan 

Proof. See Appendix E.l. □ 
Remark 8.1. FormS c > a more accurate asymptotic expression can be derived: 
m( c >~(4?V/ 2 + ^M 2 -«)/ 2 ), 



q 
p 



where 



Given Proposition 8.4, we can use Proposition 8.2 to convert the price limits into 
the volatility limits. The results are given in Proposition 8.5 below. 

Proposition 8.5. In cases (A)-(D) we have the following expressions for \i i n 

equation (8.3): 

(8.7) 





Xo-o- 


Xi 


X2 


A 


V^C^r 1 / 2 




N/A 


B 




V^C^T-1/ 2 


( 1 1 ,/^c {B A -- a '-^ 2 


C 








D 






V2n[S +< D) )r- a/2 



Proof. Follows from Propositions 8.2, 8.4. 



□ 
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(C) (C) 

Remark 8.2. For \\ , X2 more accurate asymptotic expressions can be derived: 



(C) 
Xi 



xV 



2^(^_ +C (f) T (2-«)/ 2 ) T -l/2. 



(C) 



Motivated by our discussion in Section 2, it is of interest to express these limits 
in terms of RRs and BFs. 

Proposition 8.6. In cases (A) - (D) we have the following expressions for RRs 
and BFs: 





RR 


BF 


A 




N/A 


B 




^ Cf > (-1 + 2*C<?>C<f> +4* (4f ) 2 ) r«'"V 2 


C 




16 °£ T 


D 


_ MD) _(2-a)/2 





Proof. Combining equations (2.2) and (2.3) we easily get 



(8.8) 
(8.9) 



RR(t) ~ - J -X0X1T 



-1/2 



BF(t) 



32 



X (2XoX 2 - X0X1 + 4x?) ■ 



The Proposition now follows directly from Proposition 8.5. □ 

Remark 8.3. For RR {C \ RR^ D \ BF (B \ BF {C ^ more accurate asymptotic ex- 
pressions can be derived: 

1 



RR^ 

BF (B) 
BF (C) 



-V27T 3 



2tt 
1 



/27T 



+ C ,W r (2-a)/2 \ c ,(D) T (2-a)/2 j 



27r3 C (B) 



16 

/2^3 



371 



T a'-l/2 



(Cf'+a 2 ^^ 1 -)/ 2 )^/ 2 . 



Examination of the table in Proposition 8.6 shows that in all cases (with the 
exception of the degenerate case (A)), RRs and BFs go to zero for small r. In other 
words, TSPs cannot produce finite limits for RRs and BFs for r — > 0. On the other 
hand, the rate of decay is generally slower than for regular diffusions. 

8.2.4. Tests. In order to study the quality of the above asymptotic formulas in 
detail we start with a representative LGP considered in Proposition 4.3. As in 
Section 7.3.2, we choose for illustrative purposes 1? = 0.0075, k_ = 1.00. Straight- 
forward calculation yields c£ A) = 0.0507, cffl = -0.5, = 0.0. The quality of 
the corresponding asymptotic formulas vs. exact analytical expressions calculated 
in Proposition 4.3 is shown in Table 4. Here and below C and C stand for the 



EXPONENTIAL LEVY PROCESSES 



47 



values calculated numerically and analytically, respectively, and similarly for other 
quantities of interest. 



logio ( T ) 


-2 


-4 


-6 


MA) /MA) 

^m l^m 
MA) 


0.941 
0.941 
0.304 


0.994 
0.994 
0.030 


0.999 
0.999 
0.003 



Table 4. Analytical vs. numerical values for C% , , C\ 



Next, we up the ante and consider realistic TSPs. Table 5 lists their parameters 
as well as the values of Cz, Cm, C<n as computed by Proposition 8.4. The quality 
of the corresponding asymptotic formulas is summarized in Tables 6, 7. 





a 


c+ 


c_ 


K+ 




(7 


c £ 


Cot 


Cm 


A 


0.66 


0.1305 


0.0615 


6.5022 


3.0888 


0.00 


0.1863 


0.5000 


0.0000 


B 


1.50 


0.0069 


0.0063 


1.9320 


0.4087 


0.00 


0.0670 


0.0096 


3.6492 


C 


0.66 


0.0521 


0.0245 


6.5022 


3.0888 


0.10 


0.0599 


0.1353 


-2.2867 


D 


1.50 


0.0028 


0.0025 


1.9320 


0.4087 


0.10 


0.0192 


0.0052 


-0.9610 



Table 5. Parameters for four representative TSPs and analytical 
values for C^, Cm, C<ji- Parameters for processes (A) and (B) 
are approximately calibrated to the market vols for 2y options; 
parameters for processes (C) and (D) are chosen in such a way 
that the market ATM vol for 2y options is recovered. In addition, 



D^> = -0.2395. 





log (t) 


-2 


-4 


-6 


-8 


-10 


A 




0.74 


0.94 


0.99 


1.00 


1.00 


A 


MA) /MA) 
^OT / OT 


0.60 


0.92 


0.81 


0.76 


1.28 


B 


ci B) /c (B) 


0.92 


0.98 


1.00 


1.00 


1.00 


B 


M B ) /M B ) 

°OT / U OT 


0.43 


0.86 


0.97 


0.99 


1.00 


B 


M B ) iM B ) 


1.00 


1.00 


1.00 


1.00 


1.00 


C 


cf ] ic^ } 


0.51 


0.76 


0.87 


0.95 


0.98 


C 


MC) iMC) 

^OT / OT 
OT / OT 


0.29 


0.63 


0.82 


0.92 


0.96 


C 


0.86 


0.96 


0.99 


1.00 


1.00 


C 


MC) iMC) 
^vi 1 


0.75 


0.94 


0.99 


1.00 


1.00 


D 


C^/C^ 


0.83 


0.94 


0.98 


0.99 


1.00 


D 


MD) ,MD) 
^OT / OT 


0.27 


0.72 


0.91 


0.97 


0.99 


D 


MD) /MD) 

/ ^<n 


0.91 


0.97 


0.99 


1.00 


1.00 



Table 6. Analytical vs. numerical values for Cm, Cyi, and D, 
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To numerically examine the performance of the short-time asymptotic expan- 
sions, Table 7 compares the asymptotic results given by (8.7) with the results for 
Xo — cr, Xi, and X2 obtained by a direct computation of the LL integrals in Equa- 
tions (5.1), (5.4), (5.5). Note that high-precision numerical computation of these 
integrals is a rather delicate affair; to ensure stable results, we used adaptive Gauss- 
Kronrod quadrature combined with both a judicious choice of integration region 
and a very high number of integration nodes (often exceeding 10 5 nodes). Table 
7 provides proof for the validity of the results given by (8.7), as the exact and 
asymptotic results converge to each other for sufficiently small r. Unfortunately, 
this convergence is rather slow and only truly satisfactory for r around 10 -4 years, 
i.e. for time-scales in the order of minutes or seconds. While the asymptotic ex- 
pressions do provide the correct limiting behavior, it is therefore questionable how 
useful they are for practical work. 





log (r) 





-2 


-4 


-6 


-8 


-10 


A 


log (Xo - °) 


-0.92 


-1.46 


-2.36 


-3.34 


-4.33 


-5.33 


A 


log (Xo - °) 


-0.33 


-1.33 


-2.33 


-3.33 


-4.33 


-5.33 


A 


log (Xi) 


-1.34 


0.88 


2.06 


3.01 


3.98 


5.20 


A 


log (Xi) 


0.10 


1.10 


2.10 


3.10 


4.10 


5.10 


B 


log (x ^ 0-) 


-0.91 


-1.14 


-1.45 


-1.78 


-2.11 


-2.44 


B 


log (Xo - <*) 


-0.77 


-1.11 


-1.44 


-1.77 


-2.11 


-2.44 


B 


log (Xi) 


-1.87 + TTi 


-0.98 


0.32 


1.37 


2.38 


3.38 


B 


log (Xi) 


-1.62 


-0.62 


0.38 


1.38 


2.38 


3.38 


B 


log (X2) 


0.23 


2.78 


5.16 


7.50 


9.84 


12.17 


B 


log (X 2 ) 


0.51 


2.84 


5.17 


7.51 


9.84 


12.17 


C 


log (Xo - 0-) 


-1.57 


- 2.11 


-2.94 


-3.88 


-4.85 


-5.83 


C 


log (Xo - °-) 


-0.82 


-1.82 


-2.82 


-3.82 


-4.82 


-5.82 


C 


log (Xi) 


-2.75 


-1.00 


-0.67 


-0.55 


-0.51 


-0.49 


C 


log (Xi) 


-0.47 


-0.47 


-0.47 


-0.47 


-0.47 


-0.47 


C 


log (X2) 


-0.25 


1.71 


3.00 


4.10 


5.14 


6.16 


C 


log (X2) 


1.18 


2.18 


3.18 


4.18 


5.18 


6.18 


D 


log (x - °-) 


-1.56 


-1.90 


-2.34 


-2.83 


-3.32 


-3.82 


D 


log (Xo - °-) 


-1.32 


-1.82 


-2.32 


-2.82 


-3.32 


-3.82 


D 


log (Xi) 


-2.42 + iri 


-1.95 


-1.03 


-0.43 


0.10 


0.61 


D 


log (Xi) 


-1.88 


-1.38 


-0.88 


-0.38 


0.12 


0.62 


D 


log (X2) 


-0.36 


1.63 


3.30 


4.86 


6.37 


7.88 


D 


log (X2) 


0.38 


1.88 


3.38 


4.88 


6.38 


7.88 



Table 7. Analytical and numerical values for log(xo _ c), 1°S (Xi); 



and log (% 2 ) as functions of log (1/r). 

To further illustrate the convergence of exact and asymptotic results, Figure 10 
contains log-log plots of the short-time behavior of Xo ~< J , Xi) an d Xn f° r t est case 
C (1 < a < 2, a > 0). It is evident that the straight-line behavior predicted by 
(8.7) is not realized for values of r larger than, at most, a few hours. 
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Figure 10. Comparison of the exact and asymptotic short-time 
behavior of Xo ( T ) — a i Xi ( r )> an d X2 ( T )> f° r test case C. 
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8.2.5. Options on Normal Inverse Gaussian Processes . By using (4.16), we can 
calculate short-time asymptotics for call prices both directly and via the LL formula. 



Proposition 8.7. When k ^ the asymptotic behavior of the time value of the 
) am 

d?5C NIG (v,k) _ NIG 



call option SC NIG (v, k) and its derivatives is linear in time 



c^ G (k)v + o(v), k^O, 
where v = a 2 r. 

Proof. See Appendix E.2. □ 

When k = 0, we need to analyze the asymptotic behavior of of the integrals 
[, m, n given by (8.1) with E NIG (v,u) of the form 

(8.10) E NIG (v, u) = exp (vh (h - y/u 2 + « 2 )) . 

Proposition 8.8. Consider I, m, n defined by formulas (8.1) with E NIG (v, u) given 
by (8.10) and the ATM call price C NIG (v, 0) and its derivatives with respect to log- 
strike k. Their asymptotic behavior is given by 

[(t;)~c" JG (t,,o)~- * t,ln(t,) , 



m(v)~C% IG (v,0) 
n( V )~C^ G M) 



1 

'2' 
1 



TTKV 

Proof. See Appendix E.2. □ 

As for TSPs, the rate of convergence of the corresponding integrals for NIGPs 
is slow. 

8.2.6. Options on Merton Processes. As was mentioned earlier, MPs do not belong 
to the TS class. Nevertheless, they can be analyzed by the same methods. For 
brevity we assume that a = 0. As before, we need to evaluate the integrals I, m, n 
given by (8.1) with E M (v, u) of the form 

(8.11) E M (v, u) = exp (v ( ( e 9(»+i/2)-') 2 Q(«)/2 - l) + (1 - e «) ( iu+ l 



Proposition 8.9. Consider [, m, n defined by formulas (8.1) with E M (t,u) given 
by (8.11) and the ATM call price C M (v,Q) and its derivatives with respect to log- 
strike k. Their asymptotic behavior is given by 

I (v) ~ C M (v, 0) ~ v ((1 - e") + + e*C BS ( V , -q)) , 
m(^)~Cf ( V ,0)~-isign(l-e«), 
nW~Cf,(t,,0)~^(l-e'). 

where v = At. 

Proof. See Appendix E.3. □ 
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We can obtain the result of Proposition 8.9 directly from the Merton's formula. 
To see this, let C (r, kr; a) be the normalized BS price of a call option considered 
as a function of time, time-proportional log-strike, and volatility, 

CM^)=*((-^)^)-e^(-(^ 

To the leading order we can represent the normalized price of an ATM call option 
on a MP as follows 

C M (r, 0) - C (t, -At (1 - e q ) ; 0) + Are«C (t, - (Ar (1 - e q ) + q) ; V 



We have 



Xre q C I r, - (Ar (1 - e q ) + q) ; -^= I - Are 9 ($( - + ^ ) - e~ 9 * ' 9 77 



C(r,-Ar(l-e 9 );0) = (l - e - Ar(1 - e9) ) + = Ar (1 - e 9 ) + , 

Q _ V 

kV 2 

Accordingly, 

C M (r, 0) - Ar ((1 - e«)+ + e«C BS fa, 

in agreement with our previous result (recall that, from (8.2) C = I). Expressions 
for m and n can be obtained in the same manner. 

8.3. Local Volatility Processes. Consider the general local volatility case and 
assume that an underlying is governed by the SDE of the form 

(8.12) dF(t) = a? oc (F(t))dW(t), F (0) = F , 

where af oc (F) is the so-called Normal local volatility. The corresponding log-normal 
local volatility is given by o~f oc {F) /F, for brevity we denote it as o~i oc (F). We 
consider options on the underlying driven by Brownian motion with local volatility 
described by (8.12). When F t is driven by a Brownian process with local volatility, 
we call the corresponding process a local volatility process (LVP). It is well-known 
(and had been established via asymptotic methods by [53] a long time ago) that the 
corresponding short term implied volatility Oi mp (r, k), is independent of r to the 
leading order while RRs and BFs are proportional to r 1 / 2 and r, respectively. Thus, 
market observed behavior of FX volatilities (see Table 1) cannot be explained in 
the local volatility framework. Let us assume that 

fan ^ fori 

&loc (k) = r, o lmp (r,k) = Vx n - r. 

n>0 n>0 

and express x„, RRs and BFs in terms of £„ and r. 

Proposition 8.10. In the limit of r — > the implied volatility o imp $ (k) can be 
written as 

{[k dk' = 7K dK> ) ^ 7^ 0j 
Jo ^"-r^y J Fo Tfrj^T) = b (k) . 

a loc (0) , k = 0, 

Accordingly, the coefficients XrP can be expressed in terms of £ n as follows: 

(8.i4) x<°) = e 0> x (°) = |, xf- 2 ^ !^, 
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while RRs and BFs we have the following expressions: 
(8-15) RR{T)~-M^iiT 1/2 , 



(8.16) BF (r) ~ ^ (-\toti + + l& ) r. 



Proof. Formula (8.13) is given in [5], [19], [20], [40]. Rewriting it in the form 
(8.17) 



k f k d£ 



u imp 

(k) 

JO °loc (0 

and differentiating (8.17) with respect to k three times, we obtain (8.14). Substi- 
tuting (8.14) in (8.8), (8.9), we obtain (8.15), (8.16). □ 

Remark 8.4. More accurate formulas incorporating linear dependence on t are 
also well-known (see, e.g., [5], [57], [75], [58]J: 

(8.18) C7 impil (r, k) = b {k) + rh (k) , 

where 

( b 3 a (k) , / <X ioc (0)<T toe (fc) \ , , „ 

QVPs can be used as a convenient test bed for checking the asymptotic formulas 
(8.13) and (8.18). For the quadratic volatility model with negative real roots, say, 
the terms in (8.18) become especially simple: 

ac\k 

aV* J ((HW)(W)6'W \ 
26 3 (fc) 111 I e*=q 2 fc 2 I ' ' 

2 



6i (fc) 
where 



^a 3 ((l-p) 2 + q 2 ) ((l-p) 2 + q 2 ) + 4q 2 , fc = 0. 



G (k) = arctan ( — - — \ — arctan ' 



The quality of these approximations versus the exact expression obtained Propo- 
sition 5.4 is very good, as is shown in Figure 11. However, it is clear that having 
substantial risk-reversal and straddles for short maturities is not feasible in this 
setting when parameters are of order unity. 
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o3 




(b) T=5y 



Figure 1 1 . Comparison of the exact and asymptotic expressions 
for <Timp for the calibrated QVP for different maturities. Here 
"exact" denotes the exact volatility calculated by virtue of the 
formula (5.10), while er_0, <jA are given by (8.13), (8.18). 
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8.4. Heston Stochastic Volatility Processes. We can easily analyze the as- 
ymptotic behavior of the implied volatility for the Heston model either by using 
the LL formula or by direct computation, see e.g., [71], [73], [74], [75], [77], [91], [49], 
[121], [50]. Comparison of the asymptotic and exact formulas can be performed by 
computing the relevant integrals numerically, which, by now is a well-understood 
procedure, see, e.g., [61], [110]. In view of Proposition 8.2, in order to describe the 
asymptotic behavior of the implied volatility in the Heston framework, we need to 
evaluate asymptotically the corresponding integrals I, m, n. 

Proposition 8.11. The asymptotic behavior of the integrals [,m,n defined by (8.4)is 
described as follows 

(8.19) I(t) - V^(l + A 1 r + A 2 r 2 + ...), 

OE 

m (t) - V^yr- (1 + MiT + ...) , 

Azu 

n (t) ~ — = (l + vit + u 2 t 2 + ...) , 
where the corresponding A^,/^,^ are given by (E.l). 

Proof. See Appendix E.4. □ 
Proposition 8.12. For HSVPs, we have the following expressions for \% an d 

O imp {j i &) • 

(8.20) XoW~V^(l + Xo,i T + -). 



Uq 

e 2 5 



Xi ( r ) ~ ^v^o \ l + *M T + -) 



where the corresponding \i l are given by (E.2). Accordingly, <Ji mp (r, k) can be 
written in the form 

(8.21) a lmp (r, k) ~ b (k) + rb x (k) + 

where 
(8.22) 

^ if I >■_ fx \ .< I > t . IX i j 

bo (k) = V ro o ( i + -jzr H wt^l I = v w o I 1 



4^ + 24^T = ^°l 1+ 4 + 24 



t n , , — ( , Xi,iP£k X2.i£ 2k2 \ , — { , Xi.iP 1 , X2J 2 
M*) = V^Xo,i + "4^" + "j^-J =x^o^Xo,i + — + -j r 

and Z = ek/vj a . 

Proof. See Appendix E.4. □ 

By using the well-known duality properties of the Brownian motion for r — > 
and t — > 00, it is possible to analyze the asymptotic behavior of the call price and 
the corresponding implied volatility for fixed k 7^ 0, i.e., for options which are not 
ATM. It is clear that for r — > any strike k 7^ is located "far away", so that we 
can analyze the corresponding price via asymptotic methods, more specifically via 
the saddlepoint approximation. 
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Proposition 8.13. Consider y e (— Y + ,Y_), where 
(8.23) Y± = «T*±2* 

and 4> = arctan (p/p), < 
s o (j/) - s i (y) as follows 



pe 

and 4> = arctan (p/p), < <j) < it, and define S (y) , 5i (y) = E' (y) , S i (y) 



- , x tn ysin(A'_) 
"0 (y) 



Si (y) = 



Soi (y) 



esin(X + ) 
w (sin (X_) sin + \p 2 ey) 
~e sin 2 (#+) ' 

((y- l)sin(A'_)sin(A' + ) - |p 2 ey) 



e sin 2 (X + ) 

where X± = (pey + <f> ± <f)) /2. Then for k = — Si (y), k ^ 0, the corresponding 
&imp (t, fc) can 6e written in the form 

Oimp (r, fc (y)) = (o (y) + ai (y) r + ...) 1/2 = fe (y) + &i (y) t + 
_ffere i/ie leading order terms a (y) , b a (y) are yiwen fry 

(R9d) „ u,\- fc2 (y) _ ^ 2 (y) 

& o (y) = ao /2 (y) , 

w/iere Z = ek/ru , while the higher order terms a\ (y) , b\ (y) are yiwen fry (E.5). 
Proof. See Appendix E.4. □ 
In practice, it is convenient to solve the characteristic equation 

Hi(y) + fc = 0, 

for y, and express Oi mv (r, k) in terms of k, or, equivalently, I = ek/zu . Define the 
function 3 (I', p) such that z* — 3 (Z; p) solves the equation 

(sin (XI) sin (^;) + ^) 

sin 2 (**) ~ 

where X± = (pz* + <f> ± (j>) /2. Then y* = 3 (Z; p) /e solves the characteristic equa- 
tion. While 3 does not have a simple analytical form, it is very easy to calculate it 
numerically. The corresponding calculation is particularly efficient since paramet- 
rically 3 depends only on p. We can define the following function 

, , , , ( sin (X* ) \ 

?)(M=3(M (rf^ + 'J- 

Then, for I ^ and r — » the corresponding implied volatility behaves as follows 

<Wt,0 = -t==(1 + O(t)). 
V— 22^ (t;p) 

While the above expression is not defined for I = 0, it is easy to see that 

Jwq \l\ 



V-22) (lip) 



-4- vtn , 
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so that 

T — >-U 

as expected. Moreover, it is not difficult to show that 



V-W {hp) 4 24 

in complete agreement with formulas (8.22). 

The quality of the above approximations is illustrated in Figure 12. As we 
can see, accuracy is acceptable but not particularly good. Specifically, for HSVPs 
short-time asymptotics work well for very short maturities (r <~ lm), but loose 
their accuracy for moderate maturities (r ~ 6m). 
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Delta 

— exact a_0 — o_l a_hat_0 — cr_hat_l 



(a) T=lm 




4% -I 1 , , , 1 1 , 1 , 1 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Delta 

exact a_0 --- a_l o_hat_0 — a_hat_l 



(b) T=6m 

Figure 12. Comparison of the exact and asymptotic expressions 
for <Ti m p for the calibrated HSVP for different maturities. Here 
" exact" denotes the exact volatility calculated by virtue of the LL 
formula, while cr_0, a 1 and tr 0, a \ are given by a = bo, a 1 = 
bo + rbi, and er_0 = bo, <j 1 = bo + rb±. 
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9. Wing Asymptotics 

9.1. General Remarks. It is well known that the high frequency asymptotics for 
Fourier integrals cannot be obtained in closed form unless some assumptions about 
smoothness or analyticity of the integrand are made. For instance, if / (u) is n times 
differcntiable and decays at infinity sufficiently rapidly, one can use integration by 
parts and show that 



(k)= f(u)e-* ku du = 0(\k\- n ). 



Moreover, if / (it) is meromorphic only in a strip — Y + < Imu < Y_, Y± > |, then 
it can be shown that 

(9.1) I(k) ~ n±(fc)+c±e- Y± l fe l, 

k— >±oo 

where TL± (k) represent contributions of the poles in the lower (upper) half-stip, 
respectively, and c± are constants, see, e.g., [87]. Below we show how to use this 
result in order to calculate the wing asymptotics for the implied volatility. 

9.2. Generic Exponential Levy Processes. It is easy to apply these general 
results in the case of ELPs in order to calculate the wing asymptotics of the implied 
volatility. 

Proposition 9.1. Assume that the function E (r, u) in the LL formula is analytical 
in a strip — Y + < Im (u) < Y_, Y± > |. 10 Then the asymptotics for the time value 
of a call option 8C (t, k) has the form 

(9.2) <SC(t,A;) ~ c±e-( Y±T 5) |fe| , 

k— }-±oo 

where c± are (positive) constants in (9.1). Accordingly, 



P±\k\ 



(9.3) v imp {T,k) - P±\k\, o- imp (T,k) - 

k— >±oo k— ^±00 V T 

where 

P± = 4 (Y± - y f W±)) , 0<I3 ± <2. 
Proof. See Appendix F.l. □ 

This result is well-known. For instance, Lee derived it by studying moment 
explosions, see [69]. However, its simple mathematical nature is seldom emphasized. 

We can use these formulas in Proposition 9.1 to obtain expressions for Oi mv (r, A) 
for deep OTM (A -> 0) and deep ITM (A -> 1) options. 



10 According to the Lukacs's theorem (see, e.g., [81] for a discussion), S (u) is singular at 
u = ^fiY±. 
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Proposition 9.2. Assume that the assumptions of Proposition 9.1 hold. Then, in 
the OTM case we have 



0~imp fc) ^ 
k — yog 



r 

(l-g + /2) 2 fc 



fc^oo ^ J ) k^ (1-^/2)^ ' 

mp[ ' j a^o (l-|/3 + ) 
Similarly, for ITM, we have 



Pimp (t , A) 



/3_ / 21n(l-A) 

a~i (i^i^y" 



Combination of these formulas yields 



i a\ ^ign(i/2-A) J 2 In (i - sign (1/2 - A) (i- A)) 

7 imp V- 7 " 5 



A(l-A)->0 f-t T-a \ 
2Psign(l/2-A) J 



9.2.1. Specific Exponential Levy Processes. For specific ELPs, such as TSPs, NIGPs, 
MPs, etc., the corresponding formulas can be made explicit. 

Proposition 9.3. For TSPs and NIGPs Proposition 9.1 holds provided that the 
corresponding Y± are defined as follows: 

1 



Y±=ui, /9± = 4 (w - x) , NIG. 



Y ± = K± T-, /3±=T2 + 4( K± - x /4t«±), TSP, 



Thus, for NIGPs the wing volatility is symmetric. Moreover, for NIGPs the corre- 
sponding c± in (9.2) have the form 



C± y2^|fc| 3 ^ 

Proof. See Appendices F.2, F.3 for details. □ 
The situation with MPs is somewhat different. 

Proposition 9.4. For MPs Y± = oo, so that the corresponding strip of analyticity 
coincides with the whole axis. Accordingly, tail prices decay faster than exponential. 
The asymptotic behavior of the price and implied volatility for extreme strikes is 
given by 



<5C (r,fc) ~ c±exp , a lmp {T,k) \l ~ /Tprm > 

so i/iai /or MPs t/ie raj volatility grows slower than linearly in absolute strike. 
Proof. See Appendix F.4 for details. □ 
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Additional discussion is given in [99], [15], [16] and [55], [56] among others. Fast 
decay of the call price in the wings for MPs is in agreement with general results 
presented in [1]. 

Typical cross-sections of the corresponding volatility surfaces for r = 2y axe 
shown in Figures 13, 14, respectively. It is clear the quality of the asymptotic 
approximation is far from perfect. 




(a) TSP, a=2/3 




(b) TSP, o=3/2 



Figure 13. Comparison of the exact and asymptotic expressions 
for (Ti mp (k) of the calibrated TSPs for r = 2.0 and |fc| — > oo. Here 
and below " exact" denotes the exact volatility calculated by virtue 
of the LL formula, "approx" denotes the approximate volatility 
given by (9.3), and A denotes the BS delta. 
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(a) NIGP 




(b) MP 



Figure 14. Comparison of the exact and asymptotic expressions 
for a lmp (k) of the calibrated NIGP and MP for r = 2.0 and |fc| -» 
oo. 
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9.3. Quadratic Volatility Processes. 

Proposition 9.5. For QVPs the asymptotic behavior of prices and implied volatil- 
ity for extreme strikes is given by 



5C^ v (r,k) ~ c ± e-(^)l fe l a imp (r,k) ~ J^IM 



k— >±oo ' fc— >±oo V T 

Proof. See Appendix F.5. □ 

A typical cross-section of the volatility surface for r = 2y is shown in Figure 
15(a). Once again, the quality of the asymptotic approximation is poor. 

9.4. Heston Stochastic Volatility Processes. It is clear that in order to apply 
the general formula to HSVPs, we have to determine the analyticity strip for the 
function E{t,u) given by (5.12). Analyzing the corresponding expression term 
by term, one can show that the corresponding strip is defined by the inequalities 
— Y + < Imu < Y_, where =FY± are (time-dependent) roots of the equation 

T- (T*Y±) + T+ (t»Y±) cxp (-Z ( T iY±) r) - 0, 

closest to Tl/2, respectively 

Proposition 9.6. The function E (t,u) is analytical in the strip Im«e (— Y (r) + , Y ( 
where 

Y± (r)>Y±>i 

Y± are gwen 6?/ (7.14), an d TY± (t) are the largest negative and the smallest 
positive real roots of the equation 

(py + k)+<;(y) + (- (py + k) + ? (y)) e~ eT ^ = 0, 

where <;{y) is given by (7.16). The wing asymptotics of call prices and implied 
volatilities for HSVPs can be written as follows 

6C HSV (r,k) hZ c± e-(*±^)l fe l ^(r.fc) ~ JSS, 

k—>-±oo k—>±oo V T 



w/iere [3 ± = 4 (Y± (r) - ^i? (Y± (t)JJ, < /3 ± < 2. 

Proof. Straightforward calculation. □ 

A typical cross-section of the volatility surface for r = 2j/ is shown in Figure 
15(b). This Figure shows that the quality of wing asymptotics is satisfactory but 
not perfect for HSVPs. 
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(a) QVP 




-4.0 -3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 4.0 

k 



exact approx. — Delta 



(b) HSVP 



Figure 15. Comparison of the exact and asymptotic expressions 
for <7 imp (k) of the calibrated QVP and HSVP for r = 2.0 and 
\k\ -> oo. 
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10. Conclusions 

This paper has been dedicated to the study of implied volatility asymptotics 
for a range of processes that allows for the use the Lewis-Lipton (LL) Fourier- 
integral representation of call option prices. Of key importance to us was the class 
of exponential Levy processes (ELPs), especially the tempered a-stable processes, 
but we also discussed pure diffusion processes of the local and stochastic volatility 
types. 

As we have demonstrated, the LL representation is highly conducive for asymp- 
totic work, as well-established classical methods can be used to carry out examina- 
tions of large-time, large-strike, and small-time regimes - albeit occasionally with 
some delicacies and often involving quite laborious computations. Our work is quite 
complete in establishing the formulas that characterize limit behavior for volatility 
of the processes in question. 

While some of the technical results in this paper are known, many existing results 
in the literature have been derived using a variety of (often complex) methods, and 
our paper provides a clean unification of these results along with a variety of new 
formulas. Generally speaking these results are theoretically appealing and provide 
definitive answers to a number of questions, e.g. whether non-zero limits for FX 
risk reversals and butterflies exist for the processes in question. 

In order to establish the practical relevance for the various asymptotic results, 
we have taken substantial time to undertake numerical comparisons of asymptotics 
against the exact LL solution. The difficulty of this task should not be underesti- 
mated as the integrals in the exact LL representation become challenging to handle 
numerically in the various limits. Nevertheless, once the analysis is carried out, 
it becomes clear that the performance of the various asymptotics are definitely a 
mixed bag, and overall rather disappointing. For instance, it is clear that most 
wing (i.e. large-strike) asymptotics have domains of validity that are exceedingly 
small, rendering the results mostly useless in practice. The same holds for small- 
time asymptotics of implied volatility for ELPs, where maturities generally need to 
be much less than a day for the asymptotics to be accurate. On the other hand, 
small-time asymptotics work well for diffusion processes with local and stochastic 
volatility, where it is not uncommon to find that options with maturities of several 
years can be successfully studied in the short-time limit. This discrepancy of perfor- 
mance is perhaps understandable given that ELPs (with jumps) are fundamentally 
different from diffusions, especially when observed at small time scales. 

For ELPs, what does seem to work reasonable well in many case are the long-term 
asymptotics where the domain of validity is often respectably large (say, covering 
maturities larger than one or two years). Unfortunately, this behavior is neither 
universal nor particularly robust, and it is not difficult to find examples of ELP 
configurations which, while matching the market well, result in poor long-term 
asymptotics. For instance, we have observed that for tempered a-stable processes 
the region of validity of the long-term asymptotics is sometimes dramatically re- 
duced when the parameter a is greater than one. 

In summary, while ELPs are viable candidates for describing FX market dy- 
namic, it is relatively difficult to analyze their behavior in pertinent asymptotic 
regimes. In order to accomplish this task successfully, one has to combine analyti- 
cal and numerical methods, and even then success is not guaranteed. 
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Appendix A. Fractional PIDEs 

For a G (0, 1) we can rewrite the regularized integrals in the pricing equation 
as follows 

(s) a [°° . . „-K.v\iri„\\ d V 



{e~ K ' y V (x + sy)-V (x) + (l - e~ K ° y ) V (x)) 



r(-a) Jo v v ' ' v 1+a 

=^r°i~(v(* + s y) -v { x))*L 

r(-a)J V(x) 

= e SK * x ®^ (e~ s ^ x V) - {-s) a <V, 
where V (x) — e ~ SKsX V (x). Similarly, for a g (1, 2) we have 

= T^~aj r (V {X + SV) ~ V[X) ~ SyVx {X)) e ~ KsVd (^) 

(s) a f°° . . . „ , ^ e ~ KsVd V 



/ (s (V x (x + sy) - V x (x)) + k s {V (x + sy) - V (x) - syV x (x))) 
Jo 



r(i-a) Jo 

\a — 1 roo it .-K..7/^ .7,. I _\a— 1 



= e SK * x D* (e- SK ° x V) - (-s f 1 aK a - Y V x (x) - (s) a k°V (x) , 
where V (x) = e ~ SKsX V (x). Accordingly, PIDEs (4.8), (4.9) can be written as 

V t + jV x + \o 2 (V xx -V x ) + J2 (-*)" a s e SK ° x ®« (e^V) - ]T a,K°V = 0, 

s=± a=± 

and 

Vt+ ( 7 + E aa ' aK s~ 1 J y -+r 2 - y *)+E (-«)" a.**"'*©? (e—^V)-^ a.*ZV = 0. 

\ s=± / s=± s=± 



Appendix B. The Green's Function and Call Prices for Gauss-Levy 

and Tempered Gauss-Levy Processes 

We start with the following assumptions: a g (0, 1), volatility is zero, a = 0, 
the process is maximally skewed to the left and there is no dumping, c + = 0, 
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c_ > 0, k± = 0. For brevity we omit subscripts where possible. The corresponding 
backward and forward equations have the form 

Vt+jVx+cf (V(x-y)-V(x))4L=0, V(T,x) = V T (x), 
Jo V 

G{ a) + tGW - c J™ (g<°0 (x + y)- (xj) ^ = 0, Q( a ) ( , x ) = S (x) . 

where 7 = — a = —T (—a) c = sec (an/2) i3 a > 0, or, in terms of fractional deriva- 
tives, 

V t + -yV x + aS)°(V) = 0, V (T, x) = V T (x) , 
G| Q) + 7 Gi Q) - (-1)" a^l (G (a) ) = 0, G (Q) (0, a:) = 5 (x) . 

Lemma B.l. Let (t,x) be the Green's function for a maximally negatively 
skewed stable process, such that 

P (X(t) e [x, x + dx}) = G (Q) (t, x) . 

Then, due to the scaling invariance of the corresponding process, we have 

G (a) (t,x) = -^g^ (£), 

where 

. = sec f=^, 

and g( a ) (.) is an appropriate positive function of a single variable, which satisfies 
the following ordinary integro-differential equation 

(B.l) (0 + jH 9 {a) (i + O ^ = 0. 

In particular, the Green's function associated with the LGP has the form 

(B.2) GW(t,x) = ±gW(0, 

where 

'"^ ' 3 (0 = 2V5F(-0 3/2 C< °" 

TTiis density is known as the Levy distribution. 

Proof. In order to calculate G^ we use the following ansatz 

G^(t,x) = ^ (^), 
where t = jt, and a' = 1/a. Its validity is verified below. Since 

we have to impose the following constraints on g: 

/CO 
9 (a) (0^ = L 
-00 



g (a) (Ode. 
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which, broadly speaking, play the role of initial and boundary conditions. We have 
= G { t a) (t, x) + 7 Gi Q ) (t, x) - {-l) a aV a + G^ (t, x) 

= -% N Q) (0 +9 {a) (o) - ^S^ (a) (0 • 

It is clear that i-dependence disappears and the integro-differential equation for 
g (£) becomes 

(B.4) &< a > (0 + 5 (a) (0 + (-1)" <*®l9 {a) (0 = 0. 

Here — oo < £ < 0. If so desired, we can manipulate this equation further. Specifi- 
cally, we can integrate by parts and write 

or, after integration over the interval (— oo,£], 

which is (B.l). In particular, it is easy to check directly that for the LGPs with 
a = 1/2, c + = 0, c_ = \/2?9, the corresponding f^ 1 / 2 ) is given by (B.2). The case 
of a G (1, 2) can be studied in a similar fashion. □ 

It is clear that our calculations make sense by virtue of the fact that the cor- 
responding maximally skewed stable processes are scale- invariant. Since two-sided 
stable processes possess this property as well, the above calculations can be repeated 
verbatim for such processes. It can be shown that the PDFs of standard stable laws 
have the following support: (A) the positive semi-axis if a € (0, 1) , c_ = 0; (B) 
the negative semi-axis if a e (0, 1) , c + = 0; (C) the whole axis in all other cases. 
If a € (0, 2), c + > 0, c_ > 0, then for £ — > oo we have 

5 (Q) (£K/?)~r (a+1) - 

If a € (0, 2), c + > 0, c_ > 0, then for £ — > — oo we have 

.9 (Q) (£|a,/?)H£r (a+1) . 
Thus, for all a G (0, 2), c + > 0, c_ > 0, both tails are described by power laws. If 
c_ = (c + = 0), then the left (right) tail is described by an exponential law. 

It is clear that appropriately transformed tempered stable density is not scale- 
invariant unless it is maximally skewed. Because of that, we only deal with max- 
imally skewed TSPs. 11 Thus, we assume that a G (0, 1), a = 0, c_ > 0, k_ > 0, 
and c + = 0. For brevity, we omit subscripts. In this case, backward equation (4.8) 
simplifies to the form 

r°° e~ KV du 
V t + -yV x + c ( V (x-y)-V(x))—^ = 0, 
Jo V + 

where 7 = — (a = —T (—a) (c > 0, ( = (k + 1) Q — n a . The corresponding adjoint 

equation for the Green's function G^-^ (t,x) is 

G ( r ) + 7G (., K ) _ c ^°° ^(c,*) {x + y) _ G(a , K) = Q 

11 It is clear that two-sided TSPs cannot be treated along similar lines because of the violation 
of scaling properties. 
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Equivalently, the backward and forward equations can be written in terms of frac- 
tional derivatives as follows 

V t + jV x + ae- KX D°: (e KX V) + SV = 0, V (T, x) = v (x) , 

(B.5) 

G { ^ K) + 1 G ( C> K) -(-l) a ae KX ®%_ (e^G^'^-SG^ = 0, G (q) (0, x) = 5 (x) , 

where <5= — n a a = —T {—a) n a c > 0. It turns out that for a maximally skewed TSP 
can be converted into a maximally skewed SP via a simple scaling transformation. 

Lemma B.2. Let G^ Q,K ' (t,x) be the Green's function for a maximally negatively 
skewed tempered stable process. Then 

(B.6) G<«.«> (t, x) = exp(*i+"(*-fo-*i))) g(a) (0 , 

where 



i = sec 



/an\ a (x-(i2-n)) 



t l/a 

In particular, for a = 1/2 we have 

(B.7) G (V».«) (t, x) = «P + (^-^))) g(1 /2) K) j 

w/iere 



^ v - I v v ' ■ ' ■ v . ; 

Proof. After some tedious algebra it can be show that solution of (B.5) can be 
written in the form 

G (Q ' K) (t,x) = (-(i^Je^'+^-WGW (c^^C^a;) • 
Accordingly, we can represent G^ a ' K ^ (t, x) as follows 

K a L+K(x-C,(,) (r, — Ci 

G (a ' K) (t,x) = r . g& ( 

K ' > L l/a » y L l/a 

where satisfies (B.l) subject to (B.3). In particular, for the TLGPs we have 
(B.7). ' □ 

Next, we consider prices of options for maximally skewed stable and tempered 
stable process via the Esscher transform technique (see, e.g., [54]). Due to the 
strong localization of the corresponding Green's function, we can define its Esscher 
transform as follows 

As usual, we denote the corresponding complementary cumulative function as 



/CO 
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With this notation in mind, we can write the price of a call option on a maximally 
skewed SP as follows 

C ( t) k) = £ (e* - e fc ) ^ ff <"> ( dx 



A/a 



/oo . 



k — l 

1/a J ^717^ J ■ 

In particular, for the LGP, assuming that k < t, we can represent the normalized 
price of a call option as follows 

L 2 



e 4 « 

= e t D(w,2i)-e fe D(i;,0), 

where v = 2l 2 / (t — k). In order to compute the corresponding integrals explicitly, 
we use the change of variables £ = — l/2z 2 . For maximally skewed TLGP we can 
represent the price of a call option as follows 



/'- " , x fc e t i+ K ( a; -(^- t i))- 4(( t2 -,. 1 )-x) 



e 11 r / A e Kt2?+ i« 



2^ J l>-(02-n) \ 



= e t2 D(w,2t 2 )-e fc+tl D(w,2ii), 
where v = 2l 2 / (12 — i\ — k). 

Appendix C. The Implied Volatility Expansion 
By definition, 

C fe (r,fc) = = + 9 fe a imp - e fe $ = - ^ 

\ ®imp ^ / \ G imp 

i~ ( 1 \ ± ( ^ 1 ® imp \ ( r*2- 1 o - ® imp ipk^imp) 

Cfefc (r, k) = <p\ -- 1 — ^ a fc cr imp + d fe er„„ p ^ — 

\ Vimp 2 / V 

1 fc 2 (dk&imp) ZkdkOimp J g^fj) ^ ^ ^imj 

(7jmp ^Lp ^Lp / V ^ 



^ @ imp 

'nnp 2 y \ @ imp 2 



For fc = 0we have 

(CI) C(r,0)=«(*#V *(-*>M 
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(C2) Mt,0) = ^('%^Wt)-^ * o(r) 



2 J /Vi v 7 V 2 
' ^-^IxiW + ^CCr.O)-!), 



2tt V 8 ) 2 



(C.3) 

C fefe (r, 0) = (*>M) ( %2 (r) + X! (r) - \x (r) X? (r) + ^) - $ (- 

~ vb f 1 - nr ) (* 2 (T) - (T) * ? (T) + x^y ) + Ck 0) • 



Comparison of (8.2), (C.l) - (C.3) yields 



i-4t*W) = 1(t), 



24 



Xi (r) =m(r), 



1 - ^) (x 2 (r) - Jxo (r) X? (r) + ^ ) = ft (r) , 
XoW~I(r) fl • ' 



Xi (t) ~ m(r) 1 + 



24 j 



&(r)~-^ + ^ + il(r)tf(r) + il(r)(l + ?Mj. 
Noticing that Xi (t) = Xi ( T ) rl ^ 2 > we obtain the result. 

Appendix D. Long-time Asymptotics 

D.l. Exponential Levy Processes. Due to the fact that we have chosen a para- 
metric representation of k in terms of y, so that we know that the saddlepoint is 
located on the imaginary axis at u = iy, it is easy to perform the corresponding 
saddlepoint calculation. When — \ < y < |, we can perform parallel move of the 
contour of integration to a new contour given by the equation 

Imu = y, 

without hitting singularities located at u = ±i/2. Accordingly, we can approximate 
the price of the call option as follows 

(D.l) C (r, k (y) t) = 1 + -^e^oV(v) Co (y) (l - - Cl (y)) . 

When — k + + | < y < — | or | < y < n_ + |, in the process of moving the contour 
of integration, we shall hit the points u = —i/2 and i/2, respectively. We note that 

S+ = S =k, 5_ = S (-1) = a+ (k+ - l) a + a_ (k_ + 1) Q + 7 + S = 0, 
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so that, by virtue of the Cauchy theorem, the corresponding approximate prices 

can be written as 

(D.2) 

1 - e k ^ + J^e^^co (y) (l - ±d (y)) , \ < y < ^ + \ 
i^'Mft (y) (1 - \c x (y)) , -«+ + \ < y < 



C (r,k(y) t) 



Comparison on expressions (D.l), (D.2) and (6.6) immediately yields (7.7). 
Consider y e (— Y + , Y_) and define 

(D.3) S 2 (y) = S' '(y), S 3 (y) = 4" (») , S 4 (y) = S' "' (y) , 

(D.4) 



(d.5) Co(v) = ^EL = — =L — , 

f x = L_ _ S 3 (y) A (y) 55 2 (y) _ S 4 (y) \ 

1 [V> 2~ 2 (y) v f (y) E 2 (y) f (y) + 12H| (y) 4~ 2 (y) j 

1 / (6y 2 + |) 2y5 3 (y) 55§ (y) 5 4 (y) \ 
2~ 2 (y) \(R(y)) 2 + R (y) S 2 (y) + 12S| (y) 4~ 2 (y) J ' 



(D.6) ao(y) = 2 (s+ (y) - s_ (y)) 2 , 

2 In U' 2 (y)r(y)c (y)) 

a i (y) = 

, x 2gi (y) 
a 2 (y) = — r ^ r 
r{y) 



| (r (y) ax (y) ((r (y) Or (y) - 3) (r (y) + |) - f ) + 6r (y) + 2) 

(y) «o (y) 



mfl , / s 1/2, , , , , oi (y) , , x 4a (y) a 2 (y) - a 2 (y) 
(D.7) b {y) = a ' (y) , &i (y) = 1/2 , &2 (y) = • 

2a o (y) 8a o (y) 



where 



(D.8) s ± (y) = sign (y ± £) ^-~ 0h± (y), r (y) = i?(^|£) 
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D.2. Tempered Stable Processes. For completeness we present explicit formu- 
las for Sj: 

So (y) = \° 2 (y 2 - ^ + a * ( Ks + 3 (y - 1)) - 7 (y - \ ) ■ s - 



a-1 

-7, 



Si (y) = <T 2 y + a^2 sa s (k s + s - ^ 
2 2 (y) = cr 2 + a (a - 1) ^ a s + s ^ - ^ 
S 3 (y) = a (a - 1) (a - 2) ^ sa s (k s + s (y - 
s 4 (y) = a (a - 1) (a - 2) (a - 3) ^ a s f k s + s f y - * 



s=± 



D.3. Normal Inverse Gaussian Processes. For completeness, we present ex- 
plicit expressions for E$: 



So 


(») 


= o 2 x{x- A(y)), 




(y) 


= o- 2 xyA _1 (y) , 




(y) 


= (tWA- 3 (y) , 


S3 


(y) 


= 3aWyA- 5 (y) , 


s 4 


(y) 


= 3ct 2 ^u; 2 (lu 2 + Ay 2 



where A (y) = \] up- — y 2 . 

By using (4.18), we can calculate long-time asymptotics for call prices directly. 
For k = kv and v —¥ 00, where v — ct 2 t, we have 

Lkv) = ^— f i 7 4=4^ -dx 

e-^^i (2CUKVQ 1 / 2 (y)) 

— -ay 



Q 1/2 (y) 

z kv-*vy Ki (OCjxvQ 1 ' 2 (y)) \ 

dy ) 

-5iv(y+2ojQ 1/2 (y)) 



where x = ±2xvy. Consider the following function 
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It is easy to show that on the real axis this function has a unique global minimum 
y* such that 

y* = - — , E* = S (y*) = x 2 . 
4x 

Assuming for brevity that fe < 1/2, so that this point lies outside the domain 
of integration for both integrals in the above expression, we can use the original 
Laplace method to evaluate these integrals asymptotically and get 



\ ) V irxv R ( h 



x 1 



3^2 v / 2^0 1 / 4 (i?(fc)(l + 3x 2 )+^) 



2kR (fc) QV2 x 2 R 2 (fc) WuQ 1 / 2 

where Q = Q This expression agrees with the one obtained via the sad- 

dlcpoint method. 

D.4. Merton Processes. For completeness we present explicit formulas for E^: 

5o (y) - \o 2 R (y) + A (a (y) - 1 - (1 - e") (y - 1 

S x (y) = er 2 y + A ((-g + r? 2 y) A (y) - (1 - e«)) , 

S 2 (y) = a 2 + A (t? 2 + (-q + ry 2 y) 2 ) A (y) , 

2 3 (y) = A (—q + rfy) (3r? 2 + (-g + ry 2 ?;) 2 ) A (y) , 

S 4 (y) = A (3ry 4 + 6ry 2 (-g + rfyf + (-g + 7y 2 y) 4 ) A (y) , 
where A (y) = e -?(f-s)+^ «(»). 

D.5. Heston Stochastic Volatility Processes. Substitution of u = iy in (7.13) 
yields 

k9 

So (y) = — {py + k - <; (y)) , 



fo(y) 



R( ' ' 



C (y) = \j-p 2 y 2 + 2p«y + k 2 + * . 
Differentiation of S (y) gives explicit expressions for Si (y) , S 2 (y): 

- , x «6» / -p 2 y + pk\ „ k9 (_ 2 (-p 2 y + pk)^ 

=i (y) = — p n — , =2 (y) = -t-t p + 



Let y* (fc) is a saddlepoint such that 

Si (y*) + fc = 0. 
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A simple calculation yields 

where I = p + ek/n9. After squaring this relation, we arrive at the following 
quadratic equation 

fy* 2 - 2pky* - k 2 + K _~~f 2 = 0, 
p 1 + l z 

which has two roots 

pk±Jk 2 -p 2 (k 2 - \P)Hj?+12) 

vl = y • 

It is easy to show that the relevant root has the form 

pk - sign (0 Jk 2 p> (k 2 \P) I (p2 + P) 

y* ( fc ) = j2 . 

and 

kO 1 - 

s i,± (k) = Vl {{? + pi) y* (k) + {i-p)k)±-k. 

Thus, to the leading order, <Ji mp fl (k) is given by (7.17). 

When k <~ t, the corresponding option becomes deep OTM or ITM, with one 
notable exception, which can be described as follows. Consider the implied volatility 
as a function of A, and in particular, calculate E (r, 0.5). For brevity, denote 
S (r, 0.5) by S. It is clear that the corresponding k has the form k = S 2 /2. We see 
that S can be found from the following nonlinear equation 

which can be solved numerically via a simple Newton-Raphson method. 




Appendix E. Short-time Asymptotics 
E.l. Tempered Stable Processes. If a € (0, 1), then 

dr V ' 2tt Q (u) 

Since the last integral converges and no further treatment is needed. 
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If a e (1,2), then we need to do integration by parts first (which is possible 
because k 7^ 0): 



d_ 



C(0,fc) 



ei d ( f°° e TV ( u )-' iku 



2ir dr 



Q(u) 



du 



T = 



ei d ( f°° e™(") / e 
-d 



lk:u 



2ir dr \>/-oo Q ( u ) V — 

e rv(u) 



T=0 



d 



2-Kik dr 

eJ d 
2nik dr 

k 

2nik 

k 

e 2 

2nik 



-iku d 



Q(u) 



(/: 



2u 



TV' (it) 

Q{u) ~ Q 2 («) 



T = 



OO 

oo 



t/ (u) + tv 1 (u) t; (u) _ 2^u_(w)\ T „ ( „ ) _ ife „ 

Q(u) g 2 («) 



T = 



where the last integral clearly converges. 

We are now ready to consider the ATM case. For convenience, we summarize 
our existing notation and introduce additional notation which is used below. We 
have 



v ( u ) = Yl 



s=± 



a s k s — s iu + 



f 



1 (u) = a s {—siu) a . 



s=± 



It is clear that for real u we have 

v (u) = (p + iqsign (uj) \u\ a . 

In case (A) we use change of variables v — tu and represent the integrals I, m, n 
as follows 

[ (t) = r£ (r) , m (t) = £Dl (t) , n (r) = t^OI (t) , 
where £, !JT, 91 have the form (8.6) with 

E (t, v) = exp ^T 1_a a s (t (k s - |) - + ™7 + ^ (| + <*) j 

= exp (6»i (t,w)) . 
As before, we represent £ as follows: 

£(r)=£g (r)+£f(r), 
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where 



1 f r p 1 - e 9l( - T ^ 1 



Expansion near the origin for e — > yields 

7TT [7 Q(U) J 

^ Uo " 2 + | J 2 



while the DCT yields 



so that 



(TOO f \ 



fi(r) = -(7-+*+s)^r 

Here we use the fact that 

* [Jo Q( u ) J 

which can be verified directly. For DJl (r) we obtain 

i ( f°° e 0i(r,v) 



^ Uo « 

= / sin Hv) — 

n Jo v 

= -2 si S n (7) = Cm- 

For 91 (t) we obtain 

W(r) = ^Rej^ e 01 ^^ J ~ ^Re jjf e^duj = (5 (7) = C^ . 

In case (B) we use change of variables v = t" u, where a' = 1/a, and rewrite the 
corresponding integrals as follows 

[(t) =r a '£(r), m(r) = SW(t), n (r) = t""'^ (t) , 
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where 



£7 (r, u) = cxp I ^ a s (t" (« s 

\s=± 



2 

= exp(fl Q / (t,v)) . 

As before 

£(r) = £ (r)+£f(r) 



5 \ \ a / 



where 



■Co(r) = --He( f W , du) -+ 0(e) 



1 f r°° i _ e (p+*9)l«r 
Re ■< / - dv 



7T 







*-r(-a / )Rc{(-p- lq r'} 
r(1 " a ' ) r"'cos(a' X )^4 B) >0. 



Thus, 



For SDT (r) we have 



£ (t) = (1 - a') r Q ' cos (a' X ) = C ( f\ 



n 



1 f Z" 00 e e Q '( T ^) 1 
3R(t) = IrrJ / — r — — vcfe} 

^ Uo v 
= --Im{/ cxp ((p + ig) u Q ) — 

T"" UO V 

a' 

= arg((p + zg)) 

- -— = OT ' 
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Finally, for 91 (r) we have 



91 (r) = -Re 

7T 



7T 



DO 



-Re l I e*^dv 



o 



= —Re | J cxp ((p + iq) v a ) du| 

= T -^±Re{(- {p + i q) )-"'} 

= r(1 + a V »'cos(a' X ) = (4 B) <0. 

In case (C) we use change of variables w = r 1//2 u and rewrite the corresponding 
integrals as follows 

[(t) = r 1 / 2 £(r), m(r) = SW(t), n (t) = r" 1 / 2 ^ (t) , 

where 



4 
I 

vdv 



91 (r) = -Rc<j / £ (r, u) <fe }• , 



7T 







£ (r, «) = exp fr 1 "*/ 2 £ a s (r 1 / 2 («. - |) - ««) ° + ir 1 ^ + T (2 + 5 ) - 2 

V s=± 

= cxp (0 1/2 (r,u)) . 



As before 



where 



£(r) = £g (r)+£f(r) 



^ — ^ = 



so that 



7T 7 U Z V27T 
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Likewise, for 9Jt (r) we have 



9Jt(r) ~ - z Im{ / y - y -^dv\ 

e 2— (iw / e ~ v a ~ dv 

n Jo 



Im < 


f r 




TT j 


[io 


7T l/2 


/•OO 


7T 


7o ' 



7 „-i/2 m 2"/ 2 - 1 I>/2)g - a i- a/2 

T (T T ' 



C (C) 1/2 , ^(C) l- Q /2 



Finally, for 91 (t) we get 



+ / e 2 v a dv 

2lT<J 7T Jo 

L- + 2a/2 - 1/2r W 2 +l)P J -(l + a) r l-a/2 
27TCT 7T 



= ^+d C) r 1 -"/ 2 . 



'27TCT 



In case (D) we use the same change of variables as in case (C) and obtain the 
following expressions for the relevant integrals 



1 <3 r 1 " a/2 i)(t;)-^ 



TT J V 2 TT J Q V 



2^,2 

e 2 y (u) — 

3/2-a/2 /-oo " 



-dv 

2— a 



= _^_ 2°/ 2 - 3 / 2 r(a/2-|) P j _ fa _ lvl _ a/2 

V27T TT 

For 971 (t) we get 

T l-a/2 r yoo ^ 2 ^ 

9Jt(r)~ — Im W e ~^~v(v) 

= e 2 v dv 

n Jo 

2^T(a/2)q^ 1 _ a/2 



v 



-a t 



= r (D) l-a/2 
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Finally, for (r) we get 



1 f f 00 CT 2 „ 2 
91 (r) Rc<^ / e ~ (1 + (v)) dv 



7T 



H / e 2 v "cfr; 

T JO 



2tt<7 

1 2 a / 2 - 1 /2F Ca/2 + i") » 

L_ + i {a/z+ 2 )P a -( a +i) T i- a /2 

2lT<7 7T 



'2™ OT 

E.2. Normal Inverse Gaussian Processes. First wc consider OTM case (fc > 
0). Wc have 



-NIG 



n Jk VI + X Z V Z V ' 

™ /" ( e5 - efc " 5 ) gl (fite) dx s cgr/c (fc) 
" it z 



/fc 

since the corresponding integral converges. Likewise, 



C NIG (f>j fc) = ^gL^ (,^V) rfx 

7T J fc Va; + ^ 2 v 2 v ' 
^ _uxv r°° ^ ^ c ^ /G W) 

■K J k X 



^ 2 v k 



cf/ G («, fc) = 77 _ e 2 2 2 ^ (^Vfc 2 +^ 2 ) + c k (v, 

7TV fc + >f W V 7 

7rfc 

It is clear that higher-order expansions can be calculated without too much effort 
if required. 

For the ITM case (fc < 0) we can use call-put parity and study puts rather than 
call. The corresponding results are similar to the OTM case. 
For fc = the corresponding price is given by 



^ ^ (,,o) = - ve * 2v r y- e ~ f ^ (.v^v) dx. 

tt J a Vx 2 + x 2 v 2 \ J 



Since the integrand is potentially singular at x = 0, v = 0, we have to be more 
careful than before. We split the price in two, 

C™ M) = C§ (v) + C? (v), 

where 

13 (eS — e~%) 



71" J a \X Z + X Z V Z \ ' 

Expansion of the integrand for small x, v yields 

Q(u / — ^dx = — / — dy = — In 1+ — = -uln 

0W 7T i a; 2 + ^ 2 w 2 7T i y 2 + l 2tt V ^>™' / ^ W 
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while the DCT yields 

C? (v) =0(v), 

so that 

QNIG {v0) „* v]n (l 

IT \V 

Next, we consider 



C N k IG (V,0) = ^^11 e " § Kl ( Qy / X 2 + ^ V 2\ dx = D e {v) + DT {v) t 

7T JO VX 2 + K 2 V 2 V ' 

Di (v) = f 12 Kl L^^ 2 ) dx . 

J a VX 2 + K 2 V 2 \ ' 



where 



As before, we can show that 



Finally, 



C^ G („,0) = D§ („) + D e °»~-i 



Cf/ G («, 0) = {uxv) + Cf' G («, 0) ' 



7T TTXV 

E.3. Merton Processes. We start with the following expansion of E(t,u): 



E (r, u) = exp (\r ^(V(»+i/2)-o 2 (« 2 +i/4)/2 _ i) + (1 _ e «) + \ 

= cxp [\r (1 - e«) (in + ) exp (At (V(»+i/2)-» 2 (« 2 +i/4)/2 _ 1 
= cxp (at (1 - e«) (m + 0) (l + Ar (>(»+i/2)-. 2 (« 2 +i/4)/2 _ ^ + /) _ 

Accordingly, we can represent [ (r) as follows 

1 7 1 _ e A.(l-e^ + l/2) At 7 ( e ^ + l/2)-, 2 (« 2 + l/4)/2 _ ^ 

IM = ^— / s i du-— i 



« 2 + i 



-rfu+.... 



As we know, the above integrals can be computed explicitly. The corresponding 
calculation yields 

— ((—r + KM)-'K-M)- ] 

- -A, ((1 - eT + * (i - |) + e« (. (-? - 5) - Ij - (1 - *) | + _. 

= -v((i-^r + *(M)-^(M)-a-e«)) + .. 

= At ((l-e«)++e«C BS (,,-,))+.., 
as claimed. 
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oo 



Similarly, 

,Ar(l-e«)(i U +l/2) Xt J (e^^+V^)"^ (« 2 + l/4)/2 _ ^ 

^ y « 2 +i 



m ( r) = -y — _ — i^+^y ^ 



= -I s ign(Ar(l-e < ')) + .... 



Finally, 

oo oo 
n (r) = i- | eArd-e")^!/^^ + g y ^,(i«+l/2)-^( U » + l/4)/2 _ ^ ^ + 

— oo — oo 

= «5(Ar(l-e 9 )) + .... 

E.4. Heston Stochastic Volatility Processes. We need to compute the usual 
integrals I, m, n in the short-time limit r — > 0. In this limit we have 

.4 (T .„) __-(^ + w , + 2 ,„ ( ^ w+ ^ ) ( :7 ( - 2( " )r> )) 

= -^k6t 2 Q (u) (l + ^ (pern - k)t - ^ (2 (pern - k) 2 - e 2 Q (it)) r 2 + ... 
1 — exp (-Z (it) t) 



B(t,u) = 



T- (u) + F+ (u) cxp (— Z (u) t) 

\ T ( 1 + \ (P £iu ~ k ) T + Jq, ( 2 ^ peiU ~ k ^ ~ ^) ^ 
^ (pern — k) [(peiu — k) 2 — 2e 2 Q (it)) r 3 + ...^ , 



2 
1 

"24 



C (t,m) = -^w t (l + ^-t ) C> (// ) 4 ( r. ,-/) . 



C (r, u) = — -t Q (u) y \^nj a + - k9t j (peiu - k) 

+ ^ (txj - ^ K ^ T ^j ( 2 {peiu - k) 2 - e 2 Q (u)) r 

+ i~2 n7 ° ( p£iu ~ ^ {(p £iu ~ _ 2e2< ? 1-2 + ' 

£ (TjU ) =e -^or(l+^r)o(u) e C( T ,„) 

= e -I^or(l + ^r)Q(«) + 1 T 2 Q ^ g - I roo r (l+^r) Q(u) 

x ^— tuo (peiu — k) — - [k6 (peiu — k) + wq [—2 (pern — k) 2 + e 2 Q (it)) ) t 

— ^ ^2k0 [(peiu — k) 2 — 2e 2 Q (u)) — w (peiu — k) 

x ((peiu - k) (3vu Q (u) +2k- 2peiu) + 4e 2 Q (u)) . 

It is easy to see that the relevant integrals can be calculated explicitly because the 
nontrivial part of the integrand is divisible by Q (u) and that the corresponding 
result is given by (8.19). Moreover, 
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(E.l) 



+ 192 Q ro2 (-60 (k 2 ^ 2 - fc^) + 3w 4 + 20(1 - -p 2 )fe 2 tu , 



A 2 = -(1 + ^p 2 - ^p 4 )e 4 - 20k9(3wI + 2kw - 2(1 - b -p 2 )e 2 ) 

+ Wzu 2 (10k 2 - (1 - \p 2 )e 2 )), 

Mi = — — ( -1O/S0 + 2kw - 3w 2 + 3 ( 1 - ^p 2 ] e 2 
2Aw \ \ 4 J 

v\ = g-^- ^-2 (k0 - /troo) - h7q + ^1 - ^p 2 ^ e 2 

i/ 2 = 9 (16(36k 2 6> 2 + 3n7 4 + 2e 2 Tu 2 1 + lie 4 + I2kw {w 2 + e 2 ) + Ak 2 w 2 

6144n7Q 

- 4k9(3wI + I2e 2 + lQkw )) + 8p 2 e 2 (l98n6 + zu 2 - 89e 2 - Mkm ) + 491p 4 e 4 ). 
Expressions for \i can be derived from the the general case, 

(E.2)xo,i = 2^( 6 M-^o)-(l-^ 2 ) £ 

XM = -±_(-io«tf + 2^ + 3(l-^)^ 

x 2;1 - 12( L o (- 29 ^/ £2+30 ( 1 "^ 2 ) ,ttuo 

+ 2(- 60 (l-f, 2 ). + l 9 (l-^ 2 + ^ 4 ) £ 2 

Their calculation is straightforward and is omitted for brevity. 

As many times before, we start with the LL formula (5.11) and introduce a new 
independent variable u = t~ 1 v. A simple calculation yields 

C (r, k) = 1 - — / f 2 ,i 2 \ dv > 

2?T J_ x (v 2 + jT 2 ) 

where 

v) = = 

V ' e p + ipcoth(^)' 

and &(t,v) is non-singular when r — > 0. It is clear that this integral can be 
analyzed via the saddlepoint approximation. The phase S (v, k) has the form 

so that the corresponding saddlepoint is located on the imaginary axis. Accordingly, 
we can introduce u = iy and define 

Soi (y,k) = S (y) + ky, 
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where 

i — i / \ ^0 ysin(X_) 

-0 (y) = ■ iv \ i 

and X± = (pey + <f>±(f)) /2. Expressions for f , Si (y) = E' (y), 5 2 (y) = S' ' (y), 
Cq (y) can be obtained after some tedious but straightforward algebra 



fo(y) 



Si (y) 



e 



G7 (sin (Al) sin (Af + ) + \p 2 ey) 



sin 2 (#+) 

r ('sin f > 

S 2 (y) = w p 



2 v si n(A + ) - | cos (X+) pey) 
sin 3 (#+) 



c o (y) 



fo (y) 

For k = —Ei (y), and S i (y) — S (y) — ySi (y), to the leading order we have 
(E.3) <5C(r,fc)^ Co (y)r 3 / 2 exp^ 0l(y) ' ^ 

On the other hand, as we know, for r — > 0, 

Comparison of (E.3), (E.4) shows that the corresponding <7i mp (t, fc) can be written 
in the form 

ff.mp (t, k (y)) = (ao {y) + on (y) t + ...) 1/2 = b (y) + h {y) r + 

with ao, bo given by (8.24), and a\, b\ of the form 
(E.5) 

a,iv) = *■« = ^mM' 2a ° {vy 

Appendix F. Wing Asymptotics 



F.l. Exponential Levy Processes. First, we consider OTM asymptotics, k — > 

00. In this case, in order to find the relevant asymptotics, the corresponding Fourier 
integral initially defined on the real axis has to be pushed in the lower half-plane. 
It is clear that there is a simple pole at u = —i/2 whose contribution is equal to 

1. Assuming that the corresponding shifted Levy exponent v (u) is regular in the 
strip — Y + < Imu < 0, we obtain 

(F.l) C(T,fc) ~ C +e -( Y +-5)fc. 

k— 7-oo 

Comparison of (6.5) and (F.l) yields 

k k vk 
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It is clear that lhs and rhs terms can be balanced iff 



v = /3 + , v = fi + k, a imp ~ \ < I3 + < 2, 

T fc^co V T 



(3+k 



where 



so that 



2/3 + 2 8 V 2 



/?+ = 4(y + - v / j R(Y + )). 



A similar expression is valid for ITM case with Y + replaced by Y_. In this case 
the upper limit of regularity of the Levy exponent is given by Im u < Y_ , with the 
contribution of the pole u = i/2 being e k , so that 

C(r,fc) ~ l-e k +c_e( Y - + i)'\ 



k— > — OO 



Thus 



H + IM + !ffl.( Y . + |) w . 



t) = 0_, v = |8_|fc|, ^ ~ \/ /3 =A 0</?_<2, 

fe— > — oo V T 



1 1 /?- /„ 1 

27T + 2 + ^^( Y - + 2 



/?_ = 4 (Y_ - V^(YT)) • 



F.2. Tempered Stable Processes. All we need to do in order to be able to apply 
general formulas in the case under consideration is to establish the strip in which 
the corresponding integrand is meromorphic. For TSPs it is particularly simple. 
The corresponding strip is defined by the following conditions 



■Y. ■ Tu,„ - Y . Y -- (k±T^)Y.. 



Accordingly, 



V^, P+ = -2 + 4 (k+ - yJ K 2 + - k+) k^+^, 
Y^3M, /?_ = 2 + 4(k_ - yj K 2 _ +k_) , fc^-oo. 
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F.3. Normal Inverse Gaussian Processes. For deep OTM options (k — > oo) 
we have 

„, , N Coxve^ v f°° (e§ - e fc -i ) ifi (wVi 2 + x 2 « 2 ) , 

C(i>, k) = / = -ax 

^ Jk VI 2 + x 2 v 2 

_ u)xve^ v f°° (e* vy - e k -* vv ) K x faxvQ 1 ' 2 (y)) 



■ve* v J 



k 

oo 



(y) 



dy 



— e * 2 v I y-i ( e *"V - e k~*vy\ e -2u,*vy dy 
47F 



V2^ 2j (6J _i )k v 2 y (iS+ i 

^W^-(^)* ( 1 1 ^ 



(0+i)fc , 



^ ^.^-(ui-Dfc 

27rfc 3 x 



Thus, 



Similarly, 



- - V 27r|fc| 3 *■ 



k — >■ — oo 

Accordingly, 



C(v,k) , ~__ l-e fc + ,/^ r ^^e* 2 " + ^ + ') fe . 



Q(u) 



k±=w±^, /3 ± =4(u;-x), 
so that for NIG the wing volatility is symmetric 

F.4. Merton Processes. As we know, the price of a call option on a MP is given 
where 

£ (r, u) = exp (V (~^ 2 Q («) + A )-£«(«) _ l + (l - e «) + 0^ - fc ( 

We wish to use the saddlcpoint approximation in order to evaluate the correspond- 
ing integral when |fc| — > oo. It is clear that this approximation cannot be applied 
directly since the exponent does not have a suitable form. In order to rectify the 
situation, we introduce a new independent variable v, such that 

v 

u = - — t— — — ih (k) , h (k) — > ±oo, 
h (fc) fc^±oo 



and notice that 

Q (u) ~ — 2w — h 2 , — k ( iu — - j ~ — fch, 

r (-^ 2 Q («) + A - 1 + (1 - e«) (in + 0)) ~ Are^e" 2 
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Thus, if we choose h in such a way that 

we shall be able to balance terms, and represent the integrand (to the leading order) 
in the form 

exp ^ (v) -kbj, S (v) = \ T e r ' 2lv - iv. 

We can introduce introduce y, such that v — iy, and obtain the following expression 
for S on the imaginary axis: 

S (y) = Xre-^y + y, S' (y) = -Xr^e^y + 1. 

Thus, the saddlepoint is located at 

In (Arr/ 2 ) 



and 



y = 2 



SC (t, k) <~ exp(— fch). 

|fc|— >oo 



In order to compute h, we need to solve (F.2). It can be checked that 

h ^ sign(fc) f/^m ^W^lA „ sign(fc) V21n|fc| 
|fe|~oo 77 ^ V lV 2 v / 21np-j J \k\Z00 77 

is the required solution, so that 

t-f~ 1 7 / 72h[§^\ 
oC(r, k) ~ exp — , 

Ife^oo V ^ / 

as stated. 

F.5. Quadratic Volatility Processes. In order to describe the wing asymptotics 
for the QV model we consider formula (5.10) in some detail. We need to analyze 
the asymptotics of 

c (I, K) = C sin (2kiX K ) - 2kiQ c cos {2k t X K ) , 

in the large and small strike limits. We start with deep OTM options and assume 
that K — »■ 00. For K — > 00 and l~la simple but very tedious calculation yields 

c(l,K) ~ (-l) ,+1 fc,Vp*+^q. 

As usual, large I are not relevant for our purposes since the corresponding expo- 
nential factors decay very rapidly. Substitution of this expression in formula (5.10) 
yields 

"FT" kZJ + F t (p sin (X F ) + q cos (A») ^ i?(fc ; ) sm ( 2 W 

=c+, 

This formula shows that, in general, in the large strike limit the value of the call 
option converges to a constant. For the calibrated parameters a = 1.322, p — 0.967, 
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q = 0.301, we have c + = 0.003283. Next, we consider deep ITM options and assume 
that K -> 0. For K -> and I ~ 1, we have 



if->0,Z~l 



so that 

Cft;r,g) tfq ^ kiT±*M . 

F t kZo F t (p sin (X F ) + q cos {X F )) R{h) 1 1 Fj 

if if 

=1 h c_— . 

F t Ft 

For the calibrated set of parameters we have c_ = 0.004896. These are the asymp- 
totics we need. 
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